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One of the most striking effects of the publication of Ein- 
stein 's papers on generalized relativity and of the discussions 
which arose in connection with the subsequent astronomical 
observations was to make students of physics renew their study 
of mathematics. At first they attempted to learn simply the 
technique, but soon there was a demand to understand more; 
real mathematical insight was sought. Unfortunately there 
were no books available, not even papers. 

Dr. Murnaghan's little book is a most successful attempt to 
supply what is a definite need. Every physicist can read it with 
profit. He will learn the meaning of a vector for the first time. 
He will learn methods which are available for every field of 
mathematical physics. He will see which of the processes used 
by Einstein and others are strictly mathematical and which are 
physical. Every chapter is illuminating, and the treatment of 
the subject is that of a student of mathematics and is not de- 
veloped ad hoc. The extension of surface and line integrals is 
most interesting for physicists and the discussion of the space 
relations in a four-dimensional geometry is one most needed. 
This is specially true concerning the case of point-symmetry 
which forms the basis of Einstein's formulae for gravitation as 
applied to the solar system. 

I feel personally that I owe to this book a great debt. I have 
read it with care and shall read it again. It has given me a 
definiteness of understanding which I never had before, and a 
vision of a field of knowledge which before was remote. 

JOSEPH S. AMES. 
JOHXS HOPKINS UNIVERSITY, 
June 1, 1921. 
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PREFACE. 

This monograph Is the outcome of a short course of lectures 
delivered, during the summer of 1930, to members of the graduate 
department of mathematics of The Johns Hopkins University, 
Considerations of space have made it somewhat condensed in 
form, but it is hoped that the mode of presentation is sufficiently 
novel to avoid some of the difficulties of the subject. It is our 
opinion that it is to the physicist, rather than to the mathe- 
matician, that we must look for the conquest of the secrets of 
nature and so it is to the physicist that this little book is 
addressed. The progress in both subjects during the last half 
century has been so remarkable that we cannot hope for investi- 
gators like Kelvin and Helmholtz who are equally masters of 
either. But this makes it, all the more, the pleasure and duty 
of the mathematician to adapt his powerful methods to the 
needs of the physicist and especially to explain these methods 
in a manner intelligible to any one well grounded in Algebra 
and Calculus. 

The rapid increase in the number of text books in mathematics 
has created a problem of selection. We have tried to confine 
our references to a few good treatises which should be accessible 
to every student of mathematics. 

Ch. V should be omitted on a first reading. In fact it is 
quite independent of the rest of the book and will be of interest 
mainly to students of Hydrodynamics and Theoretical Elec- 
tricity. There are several paragraphs in Ch. IV which may be 
passed over by those interested mainly in the application of the 
theory to the problems of relativity. For these we may be 
permitted to suggest, before taking up the subject matter of 
Chap. VII, a reference to an essay "The Quest of the Absolute" 



PREPACK 

in the Scientific American Monthly, March 
(UB1), and WM reprinted in the book w Relativity and Gravita- 
Mtmu&Co. (1921). It may be useful to add the well- 
nubia* of tiie Fteoch physicist, Arago "When in 



^ of the book was sent to the printer in June, 
1W1, rod it* {iahy in publication has been due to difficulties in 
ftepi^i&f tmixmft. la tbe meantime several important papers 
towing oo the Tteory of Rektivity have appeared; it will be 
gftlJdblrt to tdhr tbe reader to some significant notes by PainlevS 
itt the Comptaa Rendus of this year (1922). We are under a 
ddbt rf gnttitode to Dr. 3. & Ames for valuable advice and en- 
fiveoio interest. And, in conclusion, we must thank the officials 
of lift Johns Hopkins Press for their painstaking care in this 
rather difficult piece of printing, 

^ R D. M. 

0*UL<3H, IlL4KD. 

Jane, 1922. 
* Edited by J. Malcolm Bird. 



VECTOR ANALYSIS AND THE THEORY OP 
RELATIVITY 



INTRODUCTION 

Vector Analysis owes its origin to the German matbematiciaiis 
Mobius* and Grassmannf and their contemporary Sir William 
Hamilton, f Since its introduction it lias had a rather checkered 
career and it is only within comparatively recent times that it 
has become an integral part of any course in Theoretical Physics, 
It is well known that the subject was regarded with disfavor 
by many able physicists, among whom Sir William Thomson, 
afterwards Lord Kelvin, was probably the most prominent. 
The reason for this is, in our opinion, not hard to seek. Grass- 
mann, who undoubtedly had a much clearer conception of the 
generality and power of his methods than most of his followers, 
expounded the subject in a very abstract manner in order not 
to lose this generality. Naturally enough his writings attracted 
little attention and when, some forty years later, Heaviside and 
others were earnestly trying to popularize the method they 
swung to the other extreme and, in attempting to give an 
intuitive definition of what a vector is, failed to convey a clear 
and comprehensive idea. Roughly speaking their definition was 

* Mobius, A. P., Der barycentrische Calcul (1827). Werke, Bd. 1, Leipzig 
(1885). 

f Grassmann, H., Ausdehnungslehre (1844). Werke, Bd. 1, Leipzig (1894). 
Grassmann was particularly interested in the operations he could perform 
upon his " vectors " and not in the transformations of the components of 
these which occur when a change of "basis" or coordinate system is made. 
In this respect the point of view of his work will be found very different from 
that adopted here. 

J Hamilton, W. t Elements of Quaternions. Dublin Univ. Press (1899). 

Heaviside, 0., Electromagnetic Theory, Vol. 1, Ch. 3. London (1S93). 
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tit ** tftgtar * qwatity which, in addition to the quality 
if hffim WfcCB&ttde, has that of direction." The fault with 
k jf COOTS*, that it fails to explain just what is 
"having cfir^ctioQ*** Titat this idea requires ex- 
it dbt* wbeai we realise that the simple operation erf 
r arouBd a definite line through a definite angle 
A priori ** bis direction " in the same sense that an 
fair- Is not & vector whilst an angular velocity 
Ifc, 1itt% ^PEB, eodfess troobie arises when vectors are intro- 
dmd ia * ttramr making it' difficult to see their " direction " 
w& today BOOK* of the better text-books on the subject 
ol "jfm&ofc? vectors" sudbt as gradient, curl, etc., as 
9 they ar in any way different from other vectors. In 1901 
Bied mad Levi^Civita* published an account of their investiga- 
tions of " Tk& Absolut* Difermtial Calculus " a kind of dif~ 
fereotiatiou of vectors. This paper was written in a very con- 
densed form and did not at once attract the notice of students 
of Tlbeoretkal Physics. It was only in 1916 when Einsteinf 
cmlkd attention to the usefulness of the results in that paper 
that it received adequate recognition. However it seems to be 
the common opinion that the methods there dealt with (and 
often referred to as the " mathematics of relativity ") are 
extremely difficult. It is the purpose of this account to lessen 
this difficulty by treating several points in a more elementary 
and natural manner. For example, in an interesting introduc- 
tion to their paper, Ricci and Levi-Civita point out, as an instance 
of the power of their methods, that they can obtain easily, 
by means of their absolute differentiation, the transformation 
of Laplace's differential operator A 2 which in Cartesian co- 
ordinates takes the form 



Rica, <?., and Lo-Cwtto, T. M&kodes de Calcul dif^rentiel dbsolu. 
*OL Annalen, Bd. 54, p. 125 (1901). 
t Bin***, A., Die Qrvmdkge der atlgemeinen Rdatimt&tetheorie. Annalen 
6er Pfayak, Bd. 49, p. 159 (1916). 
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into any curvilinear coordinates whatsoever* ITiis trans- 
formation was first obtained by Jacob*,* and, while expressing 
admiration for the ingenuity of his method, they justly remark 
that it ia not perfectly satisfactory for the reason that it brings 
in ideas those of the Calculus of Variations foreign to the 
nature of the problem. Now by a method due to Beitramif it 
happens that this very transfonnatioii can be obtained by 
Vector Analysis without any knowledge of abeolute differentia- 
tion; the apparently fortuitous and happy disappearance from 
the filial result of the troublesome three index symbols of that 
part of the subject is thus explained. IB- addition we hope to 
make it dear that the methods of the " Mathematics of Rela- 
tivity " are applicable to, and necessary for, Theoretical Physics 
in general and will abide even if the Theory of Eelativity has to 
take its place with the rejected physical theories of the past. 

* Jacobi, C. (?., Werke, Bd 2, p. 191. Berlin (1882). 
f Beltrami, Ricerche di analisi applicata alia geometria. Giomale di mate- 
matiche (1864), p. 365. 



1, R**y rtokaat of physics knows the important role played 
f jft% *firf*oe *ad mime integrals in that subject. For 
ifao M&r Bia^ifcude toorfc is the line integral of the 
qpjtrtte jbn* and this will suggest a simple mode of 
tMtor, As, however, we shall wish to apply our 
III pit to gravitational spaces it is desirable at the 
&0t*fe *& deody as possible what we mean by the various 



tills tena is meant a continuous* arrangement or 
act of poists; a point being merely a group of n ordered real 
aiEB&bm. in our applications n is either 1, 2, 3, or 4 and the 
apace Is smid to be of one, two, three, or four dimensions respec- 

>rdered group of numbers we denote by x (1 \ x\ 
call the coordinates of the point they define. 
need be said for the present as to what the coordinates 
actually signify. A space defined in this way is a very abstract 
mathematical idea and to distinguish it from a more concrete 
idea of space in which, in addition to the above, we have a funda- 
mental concept called kngth, we may, where necessary, call the 
latter a metrical space and the former a non-metrical space. 
We use the symbol S to indicate our space, metrical or not, 
of n dimensions. 

SPBEADS IN S n 

It is possible to choose from the points of S n an arrangement 
or set of points such that any one point is determined by the 
value of a single variable. Thus if, instead of being perfectly 
independent, the n coordinates x (1) , -, x (n) are all functions 

Continuity is aaBumed as an aid to mathematical treatment. In certain 
modern theories preference is given to a discontinuous or discrete set of points. 

4 
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at a Jtftfl* independent variable, or parameter, 
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1,2, 



the point x is said to trace a curve or spread of one 
as i varies continuously from the value u to tyf. He points 
corresponding to the values i = ui* and tit = iii a> are called 
the end points of the curve and if they coincide, Le, if all cor- 
responding coordinates are equal the curve is said to be dot&L 
A spread of two dimensions in S* is sknflariy defined by 



where K* and tit are independent parameters. Here we have 
two degrees of freedom because we can vary the point x by 
varying either ui or ife. It is necessary, however, that the func- 
tions x*(ui, MJ) should be distinct functions of the parameters 
tii, tfc; the criterion for this being that not all the Jacobian 
determinants 



dx ( **> 



M-i,-,\ 

V* =!,-,/ 



should vanish identically. If this were to happen, we would not 
have two degrees of freedom but only one and the points would 
lie on a curve and not on a proper spread of two dimensions. 

Similarly by a spread of p dimensions in S n (p < n) we 
mean the locus of points x with p degrees of freedom; 

^2, , u P ) ($=!,, n) 



where not all the Jacobian determinants 



s p = 1, , n 



AND BELATivrry 



y. TO& we denote by V P (the corresponding 

_ being werietf) and we shall suppose all our F p 

to fct ** ffitooth "; V this we mean that all the partial deriva- 



du* \w= 1, -, 

TOb restriction is not really necessary but is 



OVKE A SPEEAD OF ONE DIMENSION I* 

aa oaxiered set of n arbitrary continuous functions 
JTli * - * JT of tlse coordinates a^, - *, a^ n) . (For brevity sake 
t liUI hereafter use the phrase " functions of position. ") 
Tim nmaerictl value assigned to the label r in the symbol Z r 
teQs wlidt one of the components X\ y , X n> which are ordered 
or imaged in this sequence, we are discussing. Now for any 
curve FI given by 

s 
n the differentials 



and then form the sum 2W 1 ) + X^dx ---- h Z n da: (n) which 
is, by definition, identically the same as 



If in each of the functions X t of position we replace the co- 
ordinates s (1) , - - , z (n) by their values on the curve FI 



2 ^* -^ becomes a function of MI, F(UI) let us say, and we may 

* Reference should be made to the classical paper by H. Poincarf, " Sur 
IB rtaidua des iategrales doubles/' Acta Math. (9), p. 321 (1887). 
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evaluate the definite integral SJf(*dd*i. Tbm h e*fled 
integral ol the ordered set of n fractions ol position (Xi, - , . 
over the carve. If, now, ire change the parameter i to a 
other parameter *i by means ol the equation KI a *fti(*i) 
points on the curve are given tyy a4& m **>(%) a f^fa) 
(^ = 1, * , n) and it is conceivable that the value of the iBtq 
might depend not only on the curve hot on the parameter 
in specifying the curve. However this is not the case since 



and 

r 



<*> 
* 



This independence, on the part of the integral, of the accidental 
parameter used in describing the curve allows us to speak of the 
integral as attached to the curve and the symbol f^&^iXjk^** 
is used since it contains no reference to the parameter n, 

In what follows we shall adopt the convention that when a 
literal label occurs twice in a term summation with respect to 
that label over the values 1, - , n is implied. Thus our line 
integral may be conveniently written 

/i = fXjtf* 

Such a label has been called by Eddington a dummy label (or 
symbol) of summation. We prefer to adopt the term *' umbral " 
used by Sylvester in a similar connection; the word signifying 
that the symbol has merely a shadow-like significance disappear- 
ing, as it does, when the implied summation is performed. 

2. INTEGRAL 7 2 OVER A SPREAD T 2 OF TWO DIMENSIONS 

Consider a set of n 2 ordered functions of position (to indicate 
which we use two labels s\, $ 2 ) 



AND 

! to *i and a* tell which one of the 
; ifiah to discuss. It is convenient to think 
if ti* fcortboft is amtaged In a square or " checkerboard " 
i n iws wad m oofemBs; then *i may indicate the row and 
Ft fe specified by means of two parameters 
i <*> (iii, *). Substitute these 
i foe the coordinates in the functions X^ and con- 
Afar tfe 4&ite dodble iirtegral 

uidu* (si and ^2 umbral labels) 

i whidi specify the points of V%. 
\ niB depend for its value not only on the spread F 2 
bet INI the parameters %, % used to specify it unless the set 
Le., X^ = X H , , t which implies the 
the n functions Xi, i; X n , n and the 
in pairs of the remaining n 2 n so that 
- z)/2 distinct functions in the set. Grouping 
^unctions of each pair we have 



T rv , , 

Is = fX^ 9 * ' dutduz (si < * 2 ) 

where now the umbral symbols do not take independently all 
values from 1 to n but only those for which the numerical value 
of ii is less than that of $2. If a change of parameters is made by 
means of the equations 



where ui and j are distinct functions of vi and v 2 the coordinates 
arc given by equations 

*<> a xW(t* lf in) s *W( fll , ^ (*=!,..., n) 
and the value of h when the b ^ are used as parameters is 
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which, by the rule for mohipiying JwxJbiaas, 



and this by the formula for the change of Variables m a double 

integral 



Starting, then, with an alternating set of functions of position 
-3Ti*i we can form an integral, (over any Ft), which depends in no 
way on the parameters chosen to specify it. To avoid all refer- 
ence to the accidental parameters we write I* in the abbreviated 
form f \X^ H d(a***>, *<*>) } (*i < ^). We adopt this in pref- 
erence to the customary notation J*{X^dxl*ddx{**>\ (*i < $$) 
since no product of differentials, such as will occur later when 
we use quadratic differential forms, is implied. 

In an exactly similar way an integral I p over a spread V p of 
p dimensions (p ^ n) is defined.* By an alternating set of 
functions X^ ^ t ..., ^ of position we mean that a single inter- 
change of two of the labels merely changes the sign of the func- 
tion. If, then, two of these labels are the same the function 
must be identically zero. Then 



is a definite multiple integral of order p extended over the values 
of u\ y - -, Up which specify the points of IV We write 



where, in the summation with respect to the umbral s}Tnbols, 
*i, ^2, , SP, *i< 82 < ' < SP- To emphasize the fact that 

* When p = n it is customary to use the phrase region of S n m preference 
to spread of n dimensions in S*. 
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f y <kw BO* <&peB$ in asy way on the parameters u ly -, u p 
fc*ffibe*rftteQ 



4 ^ x> op* = y, ar <8) 
Ji-J,J s - F,etc. 



2, 

+ Xzud(y, z, t) 



ife It we may write -T tt d(s, x) instead of JFi, adfe 2) since 
XH = Jf ij and d(z, a;) = d(x, z) 

*xample of / we may take the case of a moving 
y Euclidean space of three dimensions, the curve 
., v ~~ to change in a continuous manner as it moves, 
a:, y, 2 may be rectangular Cartesian coordinates and t 
may denote the Newtonian time. u\ is any parameter which 
serves to locate the points of the curve at any definite time 
t = t* and t may well be taken = t. Then the equations of 
our F* are 

and the parameter curves u^ = constant are the various positions 
of the moving curve, whilst the curves Ui = constant are the 
paths of definite points on the initial position of the moving curve. 
Denote dxfit by x and we have 

d(yj z)^ * 
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(It may not be superfluous to point out that it is esuential to the 
argument that Ui and t*j should be independent variables. Tbm 
in the present example t*i could not stand for the are distance 
from an end point of the moving curve if tke cum deforai* at it 
mow* although it could conveniently stand for the initial axe 
distance.) Our Is may here be written 






showing it in the form of a time integral of a certain line integral 
taken over the moving curve. Before proceeding to define the 
idea of vector quantities it is necessary to make one remark of a 
physical nature. We have written expressions of the type 

Xjlx (s an umbral symbol) 

and regarded the separate terms of these expressions Xidx, - , 
etc., as mere numbers. To actually perform the indicated sum- 
mations it is necessary, when we apply our methods to physics, 
that the separate terms in a summation should be of the same 
kind, i.e., have the same dimensions. Thus if the coordinates 
z u) . . . x (n) are a u O f the same kind the coefficients 



.... ., 



occurring in the various integrals must all have the same di- 
mensions. 

3. TRANSFORMATION OF COORDINATES 

It has already been seen that if the various line integrals 
under discussion are to have values independent of the choice of 
parameters (wi, - - -, u p ) care must be taken that the n p functions 
of position X^ t .... 8p which form the coefficients of the I p should 
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Left us now see what happens to these coefficients 
vfaa H8 change, for some reason, the coordinates z (1) , , <*> 
wft to s|afy tbe points of the Vp. The formulae of transf orma- 
n equations 



. , 2/<*>) (, - i, . . . , n ) 
tfae faecfeiooe ^ being supposed distinct so that the Jacobian of 



dtaa wt vanish identically. These equations may be regarded 
fat two ngys. First the y ( ' } may each denote the same idea as 

we have a correspondence set up 



bfcireca a point # and some, in general different, point x. 
Secondly the symbols ^ may have a meaning quite distinct 
fraa tbe symbols a^* } and then we have a correspondence 
between one set of coordinates y (t) of a point and another set 
joowlinates jr ( * } of the same point. It is the second way of 
~~dng at the matter that interests us and we speak then of a 
twmsformation of coordinates. (Prom the first point of view we 
would have a point correspondence.) Since the functions X M are 
distinct we can, in general, solve the equations* and obtain 



As an example take n = 3 and let x (l \ z< 2 >, x be rectangular 
Cartesian coordinates and (y<, y, yW) space polar coordinates 
in ordinary Euclidean space of three dimensions. 



sin y (2) COS y ( 
sin 2/ (2 > sin y ( 



yW = 



X 

.(2) 



* a. GW*rf-HHc*, Mathematical Analysis, Vol. 1 Ch 2 or 
S. B. t Advanced Calculus, ' " ' 
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In order to have a uniform transformation of coordinates so 
that to a given set of numbers y&, yP, tf there may correspond 
but one set a, a^, aP* and conversely it is frequently nma~ 
sary to restrict the range of values of one or the other set. Thus in 
the example chosen we put^ >0;0^y^^T;0^ y a> <2r. 
If now in 

Ii s J*X&x*** (9 an umbral symbol) 
we substitute 



JT, becomes Xtty, -, #*) say, and &**> = -^ <foi becomes 

mix 

ui fr an umbral symbol) 



and so Ji becomes 

9 d ^( d ^\ J ttl (r, ^ both umbral symbols) 
oy^ } ou\J 



where Y is defined by the equation 

Y r = X, (r = 1, ,n; s umbral) 

We shall from this on drop the bar notation above the X f which 
indicates that the substitution x (t) = ar ( * ) (t/ a) , , # (n) ) has been 
carried out. It will always be clear when this is supposed done. 
For an 7 2 we have 

/2 = J % X 9l8t d(x (ai \ z ( * s) ) (^i < s*) (si, $2 umbral) 

r f 3(^,^)1, , / ^ x u j c * 

= j j A, lSj 'A, ~ r f"WiQ"2 (^i < #2) by definition 

since the functions JiT^, form an alternating set. 
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k. 

Nor 



f ; -j ^ ** ** fa umbral) 

'>* M, ', * 



9&i du* dyfn* dy (r ^ dui du% 

(TI and r 2 both umbral symbols) 
BCSK if we define 

S/ ^^(*l) ^1(2) 

umbral) 



J$ takes toe f emit 

^ 

Now 



., . A , , , , 

*. . a~ 5Tfro W a mere interchange of the letters 
standing for the umbral symbols 
$1 and 82) 



m - X dxM dxM ^ since X i< * is alternating by defi- 
*" 3yfr-> fly> nitron) 

= y"ri. ri (by definition) 

Accordingly the set of functions F^. ri of position, defined as 
above, is also alternating and we may write 



Generalizing we may write / in the form 



(r\, -, T P umbral) and ^ < r 2 < < r p 
where the coeiBcients 7 ri ..... Tp form an alternating set of ' 
functions of position defined by the equations 

y y dx^ dx ( '*> 

* ..... '' " ..... '" 50 ' ' ' ~) (*i. -, p umbral symbols) 
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Accordingly , then, if an integral ovtr a oirve, or nwre generally 
a spread of dimensions p, is to h&ve a value independent of the 
coordinates the coefficients are completely determined in mry 
jgpfem of coordinate* once they are known in any particular 
system of coordinates. Tbe coefficients in a fine integral form 
as we shall see later a set of functions which " have direction " 
in Heaviside's sense and so might be called a vector. As, how- 
ever, the term vector is derived from a geometrical interpretation 
of the idea which loses to a great extent its significance when we 
apply our ideas to spaces of arbitrary metrical character the 
name has been changed and the coefficients of a line integral are 
said to form, taken as a group, a Tensor of the first rank of which 
the coefficients are the ordered components* To jjgtmgnifth 
between this definition and another of similar character this 
Tensor is said to be covariant. More generally the coefficients 
of an Ip, n p in number, are said to form a covariant tensor of 
rank p of which the separate coefficients X^, .... h are the ordered 
components. Knowing the values of the components X^ ..... ^ 
of a covariant tensor in any suitable system of coordinates x (t) 
the components in any other set y (t) are furnished by the equa- 
tions 

, , , , 

Iabels > 



Although not of such physical importance it is convenient to 
extend the idea of Tensor to an arbitrary set of functions of 
position A'* lf ..., * p which follow the same law of correspondence, 
when a transformation of coordinates is made, as the alternating 
set above. If we do this it is merely the alternating covariant 
Tensors which arise as coefficients in integrals over geometric 
figures. The reason for the correspondence between the com- 

* The term Tensor was used by Gibbs in another sense in his lectures (see 
his Vector Analysis, Chap. V, edited by Wilson, E. B.) and aJso with the same 
meaning as that given here by Voigt, W., " Die fundamentaJen Eigen- 
schaften der Krystaile," Leipzig (1898). Cf. Ch. IV, 4, infra, 
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in different systems of a Tensor in the general non- 
case would remain to be explained. 



4 IlfTOOBIJCHON OP COKTRA.VABIANT TENSORS 
10 tibfc 



Ji yX.<fe<*> EE y 7^<*> ( , umbral) 

the qmntities by which the components JT, of the covariant 
tenser ^f rank one aape multiplied have a law o f correspondence 
*fiaed by the equations 



SWlariy in the integral 



the factors Z rt , F rt which multiply the components Z r8 , 7 r<r 
respectively of the alternating covariant tensor of rank two 
have a law of correspondence given by the equations 



" ( b y definition) 

aMt^ (n, r 2 umbral symbols) 
nra (by definition) 

and so in general for an integral over a spread of p dimensions 
(P ^ n). These factors, regarded as a whole, are said to form 
a contratanani Tensor of the first, second, . . ., yth rank as the 
case may be. The sets introduced in this way are not, as in the 
case of the covariant tensors, alternating. Even though the 
correspondence between the two sets of functions of position 
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* and F n M " * may aot aria* k the *bon manner the 
set is said to form a contravariant tensor of rank p if the corre- 
spondence between the ordered components is defined by the 
equations 



The labels which serve to oHer the components are written 
o&oa* in the case of contravariant and bdow in the case erf co- 
variant Tensors. The following remark may be useful in aiding 
the beginner to remember easily the important equations defining 
the correspondence, lite umbral symbols are always attached 

to the x coordinates on the right. When the labels are 



on the left the y coordinates ^OTCI on *** I * ht 
Thus 

y,. r = -Xjifc ? (^ii i umbral) 

dy* r v dy^v 

whilst 



By an obvious and useful extension we can now introduce mixed 
Tensors partly covariant and partly contravariant in nature. 
Thus the set of n z functions of position X% TI form a mixed 
tensor of rank three, covariant of rank two and contravariant of 
rank one, if the correspondence between the two sets of ordered 
components is defined by the equations 

d?; ( '' l) 

fc> (^^^3 umbral symbols) 

Now when we recall that the x coordinates are perfectly 
arbitrary as also are the y's it becomes apparent that it must be 
possible to interchange the x and y coordinates in the equations 
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TJius, to give a concrete example, 
&aMt fee possible to derive from the n 2 equations 



wrve to define a covariant tensor of rank 2, the equations 



$2, h, fe all umbral) 



^ 1 ^ ^ fa umbral) issby the rule for 



composite differentiation and this, on account of the mutual 
independence of the x coordinates, is = unless i i = r\ in which 



case it = 



To conclude these definitions it will be sufficient to state that 
a single function of position may be regarded as a tensor of rank 
zero if its value (not its formal expression) is the same in all sets 
of coordinates. No labels are here required to order the com- 
ponents and the equation defining the correspondence is simply 

7- X 

Such a function of position is also called an invariant or absolute 
(or in the text-books on vector analysis a scalar) quantity. The 
reason for regarding this as a tensor (of either kind) of rank zero 
will become apparent from a study of the rules of operation with 
tensors. 



the formulae of transformation from rec- 
tangular Cartesian to space polar coordinates ( 3). 



THE TENSOR CONCEPT 
Here 



etc., and we obtain 



a v 

[Zj cos j/W cos y< + Z, cos y> sin y<*> - Z, sin 
^' aBteW 



[- Zi sin y< sin y + X, sin < 

the X's on the right hand side being supposed expressed in terms 
of the y's. If then we denote by R, Q, $ the resolved parts of 
the vector Xi, Xt, X 3 (the theory of the resolution of tensors 
will be dealt with later but we may anticipate here) along the 
three polar coordinate directions at any point 

}'i = R; Y, = y m e = r9; }' 3 = y m sin y* = r sin M> 
For a contravariant tensor of rank one we have 

a,.n> 
yd) = VCD "* __ i 

J ~ 



(A' (1) sin z/< 2) cos y + A' (2) sin y <2 > sin y> + A'"' cos 



J- ( A' (1) cos z/ (2) cos z/ (3) + .Y <2) cos z/ (2) sin #< 3 > - Z (3) sin y 



;a ? ,(3) 
} T (3) V<1) ?? __ i V(2) 

~~ + 



AND RELATIVITY 



Tikim *fee JPs on the right are supposed expressed in terms of the 
f*s, CaB the resolved parts of (Z<, JT<, X< 3 >) along the polar 
eeec&oate directions JB, 9, $ as before and we have 



r 

* A geoeealieaelt of which tiia is a q>ecial ease is given in Chapter IV. 



CHAPTER n 

THE AY/TORA OF TENSORS 

1, ELEMZNTAET EXILES FOB DESTVING AOT> OPERATING WITH 

TRNBOES 

(a) The Rule of Linear Chmbwa&m 
JS X2Z is a 



tensor of rank p + q( p=s J J' "^ and 
\ q = 0, 1, - / 



X%~% is another tensor of the same kind then the set of 
n?+v functions of position found by adding components of like 
order (that is with all corresponding labels, both upper and 
lower, having the same numerical values each to each) forms a 
tensor of_the same kind as X and X which is called the sum of 
X and X. By the phrase " of the same kind " we imply 
not only that X and X must have the same rank both as to 
covariant and contravariant character, but that corresponding 
components have the same dimensions. The proof of the state- 
ment is immediate for from the equations 



Gj:::5 all umbral) 

and a similar one obtained by writing a bar over Y and X we 
obtain by addition 



which is the mathematical formulation of the statement that 
X + X is a tensor of the same kind as both A" and X. 

If we multiply the equations written above, which express the 
tensor character of A r j::."5 by an invariant function of position 

21 
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(poeably a constant) m we have that mX is a tensor of the sam< 
character as X. Combining this with the previous definitbi 
of asam, repeatedly Applied if necessary, we have what is knowi 
m a fiaear combination of Tensors 



irito* &e 4 4, * are either mere numbers or scalar (invariant 
Tke sepatrate members of this linear combination mm 
kind. If , as a special case, fe is a negative numbe; 
then X+ (- X 1 ) is written X - X 



and la this way subtraction is defined. A tensor all of whos< 
components are zero is said to be the zero tensor. (It is im 
portent to notice that the property of having all the component 
ew> is an abwhde one; i.e., it is independent of the particula 
choice of coordinates in terms of which the components ar< 
expressed. This follows at once from the equations definini 
the correspondence between the ordered components in differen- 
systems of coordinates. The General Principle of Relativity 
merely says that all physical laws may be expressed each by th< 
vanishing of a certain tensor. This satisfies the necessary de 
mand that the content of a physical law must be independent o 
the coordinates used to express it mathematically. The fixing 
of the number of dimensions n as 4 rather than 3 and the inter 
pretation of the physical significance of the coordinates are th< 
difficult parts of the theory of relativity; the demand that al 
physical laws express the equality of tensors has nothing to d< 
with these and must be granted by everyone. Here we regarc 
an invariant as a tensor of zero rank.) Since the idea of a linea 
combination of tensors is reducible to a linear combination o 
the corresponding components it follows that the order of the 
separate members in a linear combination is unimportant. 

2. (6) The Rvk of Interchange of Order of Components. 

A specific example will show most briefly and clearly what ii 
meant by this rule. Consider the covariant tensor X TlT , of th< 
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second rank. The compooents have a definite order which may 
be conveniently specified by a square arrangement. 



ll 



X* 



If BOW we rearrange the if functions amongst the n* small squares 
in such a way that the rows and columns are interchanged, 
then this same interchange of rows and columns will take place 
in the square for any other coordinate system y. We denote 
the new ordered set by a bar thus 

X r , * S3 X 9 , r (r, 9 = 1, 2, -, 7l) 

Prom X r , < we obtain Y n by means_pf the equations of corre- 
spondence and we wish to show that Y T9 s Y tr where the Y n are 
obtained from the X rt by the same equations of correspondence. 
All we have done is to rearrange the order of summation on 
the right hand side of the equations of correspondence and the 
formal proof is very easy. 



= X pff -7-. . by definition (p and <r umbral) 

HlA*' ft 7/**' 



from definition of 



ss 7ar (from equations of correspondence). 

Combining this rule with rule (a) we derive some important 
results. Thus starting with X whose components are X TS we 
derive A" w r hose components are X T9 = X 9r and then the differ- 
ence X X w T hose components are X r9 AT ra = X r9 A r , r . 
This new tensor is alternating and an important example of this 
type will be given to exemplify the next rule. 
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& (*l Ths Rule qf the Simpk Product. 

Coeskler any two tensors not necessarily of the same Mad or ^ 
rack. IM us form the product of each component of the first I 
fetomdbooagpoBeEt<rftibe second and arrange the products in a $ 
definite order. The set of products will form a tensor whose 3 
rank is the sum of the ranks of the original tensors. Again it 1 
suffice to sliow I*ow the proof runs in a special example, ij 
JLet the tw^ tensors be X n and X T * and denote by the symbol \ 
Xjft tfefc product X m -X***. (Here n, r a , *i, * 2 have definite j 
iMmbal Tataes so that X**, defined in this way, is a single 
foaction out of a group of 4 obtained by giving n, r 2 , ^ ^ 
each D values frc^n 1 to n in turn.) We have to show that the 
gwq><rf a 4 functions ZJ& really form, as the notation implies, 
a tensor of rank four covariant of rank two and contravariant 
of rank two. To do this we have 



/ 

V 



by definition of 



(piy P2, cri, <r 2 umbral) 
w* 



P1PS y 

by definition of X% 
which proves the statement. 

It is quite apparent that X^is not the same as Z? so that 
the order of the factors in this kind of a product is important. 
Multiplication of tenors is not in general commutative. This 
remains true even when both the factors are of the same kind and 
rank Consider the simplest case where we have two tensors 
-ac and X bothoovariant of rank one. Then the product X-X 
is a tensor X rt = X r -X. covariant of rank two whilst the product 
X-X is a tensor X r . s X. 
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The difference X n XT* is again a covenant tensor of rank 
wo which is alternating since X n = X^ Since altercating 
ensors have a more immediate physical sygnifrmncsf tM non- 
Jteraating tensors it is natural to expect that this difference 
hould be more important than either of the direct products 
STr, or XT** It is what Grassinann called the outer product of 
he two tensers X,X in contrast to another kind of product whkt 
te calls il inner " and which we now proceed to discuss. 

u (d) Tke jRfo cf Composition or Inner Mv&pUcatia*. 

Let us first consider a simple mixed tensor of rank two Xrf* 
c which the equations of ccarespondence are 



symbols) 

f now we make r 2 = fi = r (say) and use r as an umbral symbol 
re get 



?he remarkable simplification on the right hand side is due to 
he results from composite differentiation 

agW a y <r) egco 



= if ^ 2 =t= ^i and = 1 if ^ 2 = ^i 

n this way we can form from a given tensor a tensor of fozrer 
ank (in this case an invariant). 

The proof in the general case is of the same character. 

Consider the mixed tensor ATl^JjiV"---! which is, as the 
ibels indicate, covariant of rank p + / and contravariant of 
ank p + q so that the equations of correspondence are 



, . , 

r here -r stands for -^7-; - 7-7 and so for the others. 
dz (r) dx (Tl) dx (r J 
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IH m&wA ^ymfaob of summation, *7F) on tlle ^ 



(TI Tp umbra 
dve applications of the results 




gmsustiiat 

da^o a^ r) 

Ss^'aa^^ 

= 1 if i 
so that 



(r, m, ^ all umbra 

giving the result that (JKJ^^^'J,) is a tensor, covariant 
rank J and contravariant of rank q. If q = 0, / = we have tl 
result that 



J> is an invariant (n - r p umbra 
explaining why we regard an invariant as a tensor of zero ran 
If now we have two tensors not both entirely covariant 
contravariant and take their simple product we have a mi* 
tensor to which we may apply the method here described ai 
obtain a tensor of lower rank. This is called composition 
inner multiplication of the twojensors. Thus starting wi 
X r and Z we obtain X; EE J r -Z' and then making r = * (i. 
picking the 7i diagonal elements or components of the tensor } 
of rank two)_and summing with respect to s we derive an i 
variant X.-X* which is the invariant inner product of the to 
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tensors. (To obtain an inner product the tensors most be 
of different character one covariant, the other 



Similarly from the two tensors of rank two J^ and JT^ 
first obtain the mixed tensor of rank 4 



and from this the scalar or invariant function of position 



Notice that in these cases the order of the factors is not im- 
portant the same invariant results if we change the order. 

5. (*) Converse of Rule of Composition. 

Again, for the sake of simplicity, let us explain this for a 
special case. We consider a set of n functions of position X r 
which has such a law of correspondence between components 
in different coordinate systems that for any contra variant tensor 
X T of rank one whatsoever the summation X r X T is invariant 
(r umbral). Then we shall prove that the set X T actually form, 
as the notation implies, a covariant tensor of rank one. 

We have 



y r - r s XrXM (by hypothesis) 

dar <0 . 
= Xt Y* ^ (since X T is contravariant of rank one) 

We now take as a special illustration of the tensor A' r that one, 
which, in the y system of coordinates, has all its components = 
save one which_is = 1, e.g., Y 8 = if s + r whilst } v = 1. 
This choice of X is permissible since we make the hypothesis 
that X is any tensor we wish to choose. And we have 

;MO 
Y,= X t ^p (t umbral) 

proving on assigning, in turn, to the label r the numerical values 
1, "-,71, the statement made. (It is apparent that instead of 
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* prfocdy arbitrary it is the same thing to say that 
X** dtaH b^ any ooe of the tensors which in some particular 
qygtasi nf ooocdiiMtfees ha^e each all but one of their coordinates 
m 0, the maainJBg ooe being = 1.) As another example of 
tim converse let us suppose that the n* functions X T * have such 



toMor of rack two (* umbral) where JF,* is an arbitrary covariant 
d tank two; we have to prove that the n? functions of 
/ actoaDy form, as the notation implies, a mixed tensor 
I and covariant of rank 1. 



(since X is covariant of rank 2) 
Now as our arbitrary tensor Jf let us choose that one for which 



unless both I = * and m = t 
~ 



Y, t = 1 and using ~ = 1 (r umbral) 



we obtain 

(<) <>> , , n 

((r, p umbral) 



proving the statement. The essence of the proof is that the 
multiplying tensor must be an arbitrary one. In concluding 
these remarks on the elementary rules of tensor algebra it may 
not be superfluous to remark that although, for example, the 
product X T = Xr'-Xtt is a definite tensor we do not introduce 
the idea of quotient X rt -=- XS. The reason for this is, of course, 
that there is no unique quotient; there are many tensors X 9t 
which when multiplied by a given tensor X r * in this way will 
yield a given tensor X n . In the algebra of tensors it is possible 
to have a product (inner) of two non-zero tensors equal to zero. 
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6. Applications of the Four Rides of fewer Algebra, 

The most useful applications of these niks w31 be fooad by 
returning to a consideration of the integrals which served to 
introduce us to the tensor idea. It will be remembered that & 
curve Fi is either open and has two end points as boundary or 
else is closed and has no boundaries; a spread F* of two dimen- 
sions is either open and bounded by one a? more dosed curves or 
doted and without boundaries. In general a spread Vj+i of 
p + 1 dimensions (p <L ra 1) is either open and bounded by 
one or more dosed spreads V P of j> dimensions or dbe dosed and 
without boundaries. When the spread F^_i is open tfaeare is an 
important theorem giving the value of an arbitrary integral I p 
extended over the closed boundaries V p in terms of a certain 
connected integral extended over the open V?+i bounded by V p . 
The simplest case is when p = 1 in which case an integral over 
a closed curve is shown to be equivalent to a certain integral 
extended over any surface or spread of two dimensions Vi 
bounded by the curve V\. This case was discussed by Stokes 
for ordinary space of 3 dimensions and the general theorem is 
known as " Stokes' generalized Lemma."* It will be noticed 
that the theorem is a non-metrical one as we have not yet had 
occasion to say anything about the metrical character of the 
space S n containing the spreads V p . We shall prove the theorem 
when p = 2 as this will suffice to show the details in the general 
case. 
Here the equations of the open V$ are 



and the boundaries will be specified by one or more relations on 
the parameters ?/i, u^, u s . If there are several distinct boundaries 
T 2 we may connect them by auxiliary surfaces T 2 so as to form 
one complete boundary. The parts of the / 2 over this complete 
boundary coming from the auxiliary surfaces will cancel (each 
* H. Poincar6, loc. cit. 
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tDoJitry ooeaecting surface may be replaced by two, infinites 
sid H is the integrals over these pairs o 
caace! each other in the limit as the surfaces ar 

The relation betweei 
tbft iMMMtam on the boundary may be 



*0d we intaxfaee two other functions i and tfc of %, *, u 
Mil that % % ti are distinct functions, and change over t< 
% % % as panaaerters. We shail suppose the parameters sue! 
thtt the equations giving the coordinates x are uniform hot] 
Not oedy cioes an assigned set of parameters give ; 

x but to a point x there corresponds but one set o 
ptnuaeters*. 

Accordingly the surfaces i = const cannot intersect each othe 
and they form a set of dostd level surfaces filling up the initia 
opa F t . On each of these closed level surfaces we shall hav< 
the level curves i = const., tfc = const., and we suppose th< 

* of tii, u*, w s so chosen that these level curve 



Now consider the integral 

It s fX^d(x^, x ( '*>) (si, 2 umbral and Si < s 2 

extended over the boundary t> 3 = 0. If, instead of integrating 
over i = 0, we take it over any of the level surfaces 2)3 = constan 
it will take on different values depending on this constant anc 
to indicate this we write 

7 2 (r 3 ) SEE 



r x ~- K 4. 

** dv* 
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(It is only necessary to differentiate the integrand since the 
limits of the integral are independent of s). Now if F is any 
function of position (not merely of the parameters) 4 ' 09 a doeed 

3F 
curve with parameter t> the integral S--d9 taken round the 



curve is necessarily zero. For it is the difference of the 
values of F at the coincident end points erf the curve. If, in 
particular, we take as F the function 

* f , ^ 

fa, * umfaral) 



and integrate round the dosed curve t = constant we get 

and integrating this with respect to ^ over the surface t> 3 = con- 
stant we have 



.dX^dafi'tdiP*} 

+ r~^-^ ^ r 

ovi dvz ofy J 
Similarly on taking 



a r 

and integrating J* -- dvidvz over the closed surface z?s = const. 
di^ 

* The distinction implied here should be clearly grasped. If the equations 
of the curve are 

Xi = a cos v 
x* = a sin v 

F must be periodic in v with period ITT. 



S3 VBOTOR ANiJraSIB AND EELATTVTTY 

wept 



a* 



Nkw add t$Mm two equations together and note that 



i in the smnmation cancel out in pairs owing to 
tbe tteroatiQg ch&reier of X,^ the factor multiplying X^ 
ia the swmmation bong obviously unaltered by an interchange 
of the symbols *i and %, We find that 



so that 
'-./ 



Now the J^ are functions of position, i.e., of the coordinates x 
so that 



The second term in rf/j/rfr, we shall slightly modify by a change 
in the umbral symbols. Thus 



= . 

at, a^ ait " d?"> "aTT ~fa~ ~ (Sl> Sil Sa a11 umbral) 
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so that we can write 

dl l= r \ aZ^ _ BX^ _ dX^ 1 3sW daf*d** ^ ^ 
^ ->/ l&e^ az<> da** J d* d* d* *"*** 

On writing 



T = d^i<* _ d3T<m _ &%** 
**** "" &e<*> dy?-** fe<*> 

and integrating the expression for dl^de* with respect to 
we find 



-T 

C/Uj 



since the set of functions -ST^a^ defined as above is obviously 
alternating (on account of the fact that X T9 is an alternating set). 
The limits for t> 3 are 7? 8 = and v* = some constant for which 
J 2 = since the corresponding F 2 is either a point or a spread 
traced twice on opposite sides. Let the integration be such 
that 7)3 = is the upper limit and we have 

7 2 = fXt lH d(x^x (9 *} (*iO 2 ) over boundary 

= fX^d(x^x^x^ (si < * 2 < * 3 ) overtheF 

In general from 

d&** - - *<>>) (*i < 8, 



j.* 



over a closed boundary we derive as equivalent to I p an 

/H.! 
where 



* It will be observed that placing the 4- sign before I t on the left makes 
vi = the upper bound of the integral f ~ dv t . Thus vi is increasing away 
from the open spread Vs. 
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ft ll usual to piutfev^ * cyclic arrangement of suffixes for the 
X* *i*i thea, OB *<xoi^ 



tbs oppr aig8 being used when p is even and the lower when p 
m ocli Sfaace J, is by hypothesis invariant so is Jp+i because 
IjM* "* ^ ^ accordingly the coefficients X^ ... Vl form an 
dfagMtiiH oovariant tensor of rank p + 1 [seen either directly 
m wben tensors were introduced or as a case of the converse of 

),th act of functions- ----- ~- - <foi - ^p+i fonn- 

0*1 



log in irbitrtry eontravariant tensor of rank p + 1]. In this 
w*y we can derive from any alternating covariant tensor, by a 
species erf differentiation, a covariant tensor of higher rank. 



p = 1. Pram any covariant tensor X r of rank one we derive 
an alternating covariant tensor of rank two 



It is the negative of this tensor that is called the curl of the 
vector A" in the earlier vector analysis. It is rather important 
to notice that this, and the other tensors of this paragraph, have 
no reference to the metrical character of the fundamental space 
S. The derivation of them by the methods of the Absolute 
Differential Calculus introduces, therefore, extraneous and un- 
necessary ideas. 

p = 2. From an alternating covariant tensor of rank two 
X n we derive the alternating covariant tensor of rank three 

V a* d-<V" , dX tt , dX tr 

r "~~ i ~ + 
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If n = 3 there is only one such function and in tiv 
it is called the divergence of X n . We shall hav 
slightly for the general tensor analysis. It is 
notice that if we take as X r9 the tensor of the j 



we find X n t = 0. It is easily seen that this happens in general. 
If we derive X^ ... Sf from X^ ... ^ in this way then the 
X* v* d 63 *^ from X^ ... h is s 0. Wfcen tiie JT % , ^ 
derived from X^ ... is = we have that I^i ss and so I 9 
(extended, of course, over any dosed spread of p dimensions) 
is s 0. In this case I P is said to be the integral of an exact 
diferential. It can then be proved that the value of I p over 
any open V p is equal to the value of a certain integral J,_i over 
the closed boundary of this V 9 .* 
*If 



is the integral of an exact differential we have f * ^J m ^ -^ -^ 

partial differential equations 

^ia p+i s ^ 
The theorem stated is that these are the necessary and sufficient conditions 

that there exist ( ^ j) functions of position X tl ... 9p _ l satisfying the ( 
partial differential equations 



That the conditions are necessary is an immediate result of a direct substitution 
of the left hand side of the equation just written for X 9l ... , p in the equation 
of definition 



To prove the sufficiency an appeal is made to the principle of mathematical 
induction. Let us, for definiteness, take p - 2. Then we shall prove the 
statement that if the theorem is true for a particular value of n it is true 
for the next greater integer value n -f 1. Granting this, for the moment, we 
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!, fife ie the BCT* and last case if n = 4. For 
of it is second in importance only to the 
"* 1* In order to avoid having to write out separately 

AM* $* tfcwwa is tree for 2. (In this case there are no ra- 
on acocxmt of the alternating character of 



varoflfamg expresses these conditions, it is neces- 



S55 



nor Zi involves 

tiflft JT! My b aay Emotion of ar and Z, - - y^-^ZuefcW, the 
(f ?JM<KH* s,^. In the integration g<*> is regarded as a 
by the iadoction kmma the theorem is true f or n = 3 
to* - 4 tad ao lor every integern. 
l^pro^ tbomdaofeioii lemma iet us seek for a solution of the equations 

X 9 . 

Jw (r 

0. We have then 



whew ^<> is a constant; X r is any function of s&>, , <-^) and in the 
istogFatkm a^, - - , x**- 1 ) are constants. The remaining equations 

cr _ 

Xrt 
grve on substituting these values 



_ ^ 
X r . - JCr. + g^ - ~ where X ff is the function X r , when z) is 



put ,<). 
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cases corresponding to n even and n odd we shall adopt the 
first form for X^ ... ^ 



Heooe we have the (* ~ l ) equations 

= a? f al, 
* "* 6x<*> ax< r> 
Kith n - 1 unknowns X r and involving n - 1 independent varabfee 

so, , <I *^ > . Also we have r w ~ J int^rabfl^y eqoatioos ZT 

found by putting s c > - *o <ll> ia 

Z r ,i = (r < < * - I, -, - 1) 

Hence if we can solve these equations, i.e., if our hypothesis is true fcr - 1, 
^e can solve the original equations which are identical in form but involve 
one more independent variable a?<*>. The particular case of this theorem 
corresponding ton = 4, p = 2, tellsus that Maxwell's equations 



l a 

curl 2?-h--rr = divB=0 (in the usual notation) 
c ot 



imply the existence of the electromagnetic potential (A x , A A^ 

which is as in the general case when p * 2 a covaiiant tensor of rank one 

such that 

B - curl 3; E = - grad $ ---jjj- 

For further details cf . Physical Review, N. S., Vol. 17, p. 83 (1921). 

It is apparent that there is a great degree of arbitrariness allowed in the 
determination of the functions X^ ... , p-1 ; in fact we may add to any solution 
any alternating covariant tensor of rank p 1 whose integral over any closed 
spread Vp-i of p 1 dimensions is zero. For example we may add to the 
electromagnetic potential any gradient of a function of position; that is 
if (A xy A tf , A,, - c<t>) is any determination of the electromagnetic potential, 
so is 



N where F is an arbitrary function of x, y, z, t 
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f Hh I a and there & only one distinct function X^ ... ^ 
itei^^ Let us choose 

formula is 



l, 2 



Now tivew re 00)7 n distinct functions Z^ ... ^ and it will be 
m& conv^siient, to indicate these by means of a single 



(X*) 



where we trc careful to put parentheses round the symbols (X T ) 
to indicate that they are not the components of a covariant 
team of rank one. 

MJUW^ avaiWhinwelf of this arttrariness and chose Fs^ s 



yioidiog, from the theory of the Newtonian Potential, 
- 1 .divA 



div A + I dj 
c at 



The usual procedure with modern writers is to choose F so that 
divl 

Ite equation determining F is now 
whence 



f(di v ^I) 

, I V ca/^r 

1 / i-dr 

4r*^ r 



from the theory of the retarded potential. 
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Then we have 

^ -Hr (* an mnteal label) 

Although the (Z ,) do not form a covariant tensor of rank one 
tb^ are very dosely related to a m^^ 
In fact there is a reciprocal relationship between an aflama&t? 
oowriani tensor of any rank r and an allied cowbraaanant alter- 
nating tensor of rank n r. It is a special case of thkredprodty 
stressed so much by Grassmann in his Ausdehmtngslehre tlmt 
gives tte dual relationship of point and plane, line and line in 
analytic projective geometry and it is from the terminology of 
that subject that the terms " covariant " and " contravariant '* 
are taken. In order to bring out this reciprocal relationship in 
the dearest manner we must make a digression and discuss what 
are meant by " metrical properties " of space. 



CHAPTER ffl 

09 TOR MBTRICJLL IDEA INTO OTTB GEOMETRY* 

Let w wonder a curve F* speriied by the equations 



*dr^di3eren^ 

$&*&** (r, * umbral) 

wfcrro the f **e functions of position, will be invariant provided 
thai these fuixrtions form a covariant tensor of rank 2. (This 
is * coosequenoe of our rule (d), Ch. 2, 4, and its converse since 
the set of a* functions 

,, , . <k< r > <c<*> , , , 2 
dsf*dx* s -- - (dw) 2 
dw du 

form a contravariant tensor of rank two.) Accordingly the g r9 
being of this kind the integral 



9 == 




du du 

has a value independent of the choice of coordinates x; it is called 
the length of the curve V\ from the point specified by u$ to that 
specified by u'. If the upper limit u' is regarded as variable 
and written, therefore, without the prime S is a function of this 
upper limit u and its differential is given by 



(r, * umbral) 

where the positive radical is taken on extracting the square root. 
It will be convenient to agree that, in some particular set of co- 
ordinates a% we arrange matters so that g r , = g. r \ this can always 

* The moet satisfactory presentation of the general idea of a metrical space 
is that given in Bianehi, L., Lezioni di Geometria Differenziale, Vol. 1, 152. 

40 
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be done by rewriting any two teams* &&&&*$ + 
for example, of the summation which do not satisfy this reqrare- 
inent in the form ifot + ^<fe<fa + \($* + &&&*&? 
The equations defining the covariant correspondence 



where 

then show that 



We may express this result by saying that the property of any 
special tensor of being symmetric is an absolute one just as is 
the property of being alternating. 

2. RECIPROCAL FORM FOB (ds) 2 
Consider the n linear differential forms 

r ss grJxW (s umbral; r = 1, , n) 

We can solve these for the differentials dx ( * ) in terms of the n 
quantities r as follows. (Note that the r form, as the notation 
indicates, a covariant tensor of rank 1 from our rule (d) of com- 
position or inner multiplication.) Let us denote the cofactor of 
any element g rs in the expansion of the determinant 



nl ' ' * ffnn 

by (G rt ), observing in passing that (<?) ss ((?). The parenthe- 
ses indicate that the (G>.) do not form a tensor. From the 
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s 35 # when m * (r umbnJ) 

when m 4= * 

We iiiiftP aw introduce the hypothesis that our metrical space 
iv m& &*t f does not vaaish i(teDdi(By ^t will be presently seen 
tfct th is *& afceoiirte property) and for all points where g is 
aotserowehave 



- (rumbral) 

when m 4 s s 

Write I?"* * {OiwJ/f &d fet us justify the notation by showing 
tb*t tbe f 1 * form a contravariant tensor of rank two. From our 
nyrintt jt is symmetrical and so we have in addition to 



=1 if TO = * 

= ifm=M 
the equivalent equations 

giro** =1 if m = s 

= if m 4 s ^ 

These relations suggest that we multiply the equations of defini- 
tion 

r . gjkM 

by f and use r as an umbral symbol. We obtain then 



(r, 5 umbral) 
== ir( w ) from our relations just written 
Accordingly 

(&y^g*iaP>daf>*> = g*^ &{?*'& (I, m, r, 9 umbral^ 
s ^ r *Sr^ (r, 5 umbral) 

since j^ s unless m = r when it s 1. 

The r. form, by rule (c), Ch. 2, 3, an arbitrary contravariant 

tensor of rank 2 and (d*) 2 being, by hypothesis, invariant, the 
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converse of rule (d), Ch. 2, 5, gives us the result that the jf 
form a contravariant (symmetrical) tensor of rank 2. When we 
write 

(&)*/&& (f,*umbral) 

it is said to be written in the reciprocal form. We could start with 
this form and write 

dx> as f*& 

and solving these obtain 



and then find 



3. If now we have two determinants a s [a), i s (&] 
each of order n (the notation implying that a r * is the element in 
the rth row and $th column of the determinant a) it is well 
known that the product of the determinants a and 6 may be 
written as a determinant c* of which the elements c r , are defined 
by 

c rt = airbis (I an umbral symbol) 



This kind of a product is said to be taken by multiplying columns 
of a into columns of 6. 

We can, with the aid of this rule, easily see how the determinant 
g behaves when we change our coordinates x to some other 
suitable coordinates y. We get a determinant / of which the 
r , 5th element is 



Here -^ may be conveniently denoted by (ji r ) since it is the 

Z, rth element of the Jacobian determinant 7 of the transformation 
from x to y coordinates 

* Cf. Bdcher, M.> Introduction to Higher Algebra, Chap. 2, Macmillan (1915). 



jUfALYHB AND BBLA-Tiviii 
3as> 



a* 



T x , 

is the. & deraaat of the product gJ so that ^w 



fe tibe rth element of ti product of the determinants gJ by J. 



important formula shows us that if g + neither will 
arieaa J s in which case the y's would not be suitable 
/ can be zero at points where g > * if J = at 



. 

thoee points; such points would be singular points of the system 
of coordinates and the quantities f* would not be defined for 

them. 

EXAMPLE 

In Euclidean space of 3 dimensions with rectangular Cartesian 
coordinates af x w a^ we write 



so that ju = ff = ffw = 
coordinates we find 



= y = - In s P ace P olar 



,(2) 



Here g - 1 

so that 

9 
1 



II-/U-/M-O. 

/ = /11/S2/33 



= 13 = ff 53 = 







sn 
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and 



IE fact i = dxP*, etc. There are DO singular points in the x 
coordinates but there are in the y system; those for which J = 0, 

ie., 

yW*siay<& = 

These are the points on the polar axis 

ytt> ss r 0; 2^> 9 = or T 

4. If now ui tt are any independent parameters in terms 
of which it is convenient to specify both the x and y coordinates 
we have, by definition of the symbol, 



and a similar equation for d(x (1) ar (n) ) so that 



If we multiply the determinants -~r ~ and (I ^ (j j 

together and note that 

j~: =1 if t = r (s umbral) 

= if t*r 

we find that their product is unity and so we can write the 
quotient 



, 

*P>) aum x x 

(m umbral) 
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so th*t this expression is an wwwrwm*. In view of the fact that 
ft depends oa the fundamental quadratic differential form (eb) 2 
it is ctfldd a mtiriexiJ iavariaat. 

Let m consider an integral over a region of the fundamental 
&* /Zi^^a^ *). Here XL.., is the ain^fe 
function of an arbitrary alternating co variant tensor of 
rak m Since the integrand is invariant and since V</ d(x (1) - * 
is invariant it follows by division that X\ ... tt -s- V<7 is an 
As an application of Stokes' Lemma we have already 
E 



(where X^...^ is any alternating covariant tensor of rank 
n 1) then 

^ umbral) 



b the coefficient of an integral over a region of S n . We see 

1 d 
therefore that-p ^ (X.) is an invariant. 

ijg dx w 

We shall now investigate the nature of the n functions (JQ. 
Under a transformation of coordinates from x to y we find, for 
example, 



(^i - s n _i umbral) 
X d(x ( '*> . . a ( *n-i)) 

<^ ^""-^(yW -.. y<-) 
(owing to alternating character of X^ ... ^ 
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And, accordingly, if we denote the cof actor of 
sion of / by (Jr.) we have 



In general 

(F r ) (/)(*) 

If we solve the n equations 



ifa + i" r = 1 



so that we may write 






(Jf ) 

showing that ^-~ is a contravariant tensor of rank one. We 

10 
may then put (X 8 ) = ^g X* and our previous result takes the 

form that -7= ( V^ X 8 ) is an invariant; -3T* being any cx>ntra- 



variant tensor of rank one. This metrical invariant is known 
as the divergence of the contravariant tensor. 

5. SPECIAL RESULTS 

If u(x w - x (n) ) is any invariant function of position the rule 
of differentiation 
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(fat 

fli i ttatt the it fuactioBs X t m -jfom a covariant tensor 

rf rmnk one; this is known as the tower gradient. If Z r is any 
wvamnt tensor of raok one its simple product by itself or 
* ttpttfe ** is * covariant tensor of rank two, X n ^ XrZ,. 
Hesccbynife(cf), Ch. 2, 4, 

fXrX* is an invariant (r, a umbral) 

Hii is caflesd tbe square of the magnitude of the tensor. In 
pttfadbar the square ol tbe tensor gradient is the invariant 

A du du , . _ % 

Al * s f a^d^y (r > * umbral) 

Has is known as the " jfr^ differential parameter of u." Similarly 
the magnitude of the square of a contravariant tensor of rank 1 
is the bmriant g^X^X^. 
Again 

g"j = X* (r umbral) 

is contravariant of rank one (rule (d)). Hence 

Id/,- du \. 
Tgtofi \ V * f9 flP J B an mvariant ^ * u 



by the result of the preceding paragraph. It is written A z u 
and is known as the "second differential parameter."* In 
ordinary space of three dimensions in which the ar's are rec- 
tangular Cartesian coordinates 

g r = if r =|= $ 
= 1 if r= * 

and g n = /; V = 1 so that A 2 w takes the form 

2 , * 

9 ~l~ 



r, J., Transactions Cambridge Phil. Soc., Vol. 14, p. 121 (1885), 
obtains this transformation in the case n - 3 by the application of the Calculus 
of Variations, 
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When we change over to any "cm 
we have under the form 



the expression of this magnitude in a io~ 
coordinates. 

6. GENERAL ORTHOGONAL 

"Whenever we have, in any space, eo 
expression (d*) 2 involves only square te 
the coordinates are said to be orthog __ 
explained later). It is usual to write, in this eat*;, 



* 
accordingly 



so that 



g= ...g m=n . f . = o f 
The square of the gradient is 






whilst the quantity 



L d / AI , 

An ] T-7TT ( T - T Tm I + 

) \ A 2 



The reader should write out the explicit formulae for space polar 
and cylindrical coordinates in ordinary space of three dimensions. 
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7* T!m SPECIAL OB KESTSICTTED VECTOB ANALYSIS 

In the form given to the theory by Heaviside and others only 
tboee oootdiMtes^ or ^ were consider^ which the fundamental 
metrical form is 



coordinates we call rectangular or orthogonal Cartesian 
coordinates and the space we call Euclidean. It is true that 
coe was made of Stokes* Lemma to find expressions for important 
fawmats as * in other than orthogonal Cartesian coordinates 
bet oo atfcesspt was made to define the components of a vector 
m these coordinates. Now when we restrict ourselves to that 
subgroup (of ail the continuous transformations) which carries us 
bom oae set of orthogonal Cartesian coordinates to another 
the distinction between covariant and contravariant tensors com- 
pletely disappears* The transformations are necessarily of the 
linear type 

(s umbral, r = 1 - n) 



where the a*s are constants. Since here / = g = 1, J 2 = 1* 
and so the equations just written have a unique solution for the 
jf's. To get this most conveniently note that dx (r) = (a r9 )dy^ 
and squaring and adding we have 

(a r ,)(a r ) = t =j= s (r umbral) 

= 1 t = s 

Hence multiplying the equations for x by a rt and using t 
as an umbral symbol we find 

(a rt )z< r) B (a rt )(a r .)2/ ( ' ) (r, * umbral) 



Accordingly the equations of correspondence defining covariant 
* We shall consider only direct transformations; those for which J = + 1. 



THE METRICAL CONCEPT 51 

and contra variant tensors are, for this restricted set of trans- 
formations, identical. Again denoting by (A n ) the cofaetor 
of (an) in the expansion of the determinant J we have by the 
usual method that 



and since the solution is unique we must have (on) 
Hence since g = 1 we have that the n distinct compaamb of an 
alternating tensor of rank n 1 form a fcnaor # raai 011*. It 
b for this reason that when n = 3 it was found necessary to 
discuss but one kind of tensor that of the first rank which was 
called a vector.^ Still some writers felt a distinction between the 
two kinds; that of the first rank they called polar and the 
alternating tensor of the second kind, whose three distinct com- 
ponents form a tensor of the first kind, they called axial. Thus 
a velocity or gradient are polar vectors (the first being properly 
contravariant, the latter covariant) whilst a curl or a vector 
product are axial vectors. 

When, in the mathematical discussion of the Special Rela- 
tivity Theory, it was found convenient to make n = 4 [the trans- 
formations (Lorentz) being still those of the linear orthogonal 
type], a new kind of tensor or vector is introduced. Here it is 
the alternating tensor of the third rank which, when we consider 
merely its four distinct components, is equivalent, from its 
definition and the properties of the transformation, to a tensor 
of the first rank or " four-vector." There remains the alter- 
nating tensor of the second rank and the six distinct components 
of this were known, for want of a better name, as a six-vector. 
As an example of the general theory we have that 

dX (t} . 

(a) the divergence of a four-vector " is an invariant. 

ox 

(s umbral) 

* This is merely a special case of the previous result that J ^ = (J tr ). 
t Until a consideration of non-alternating tensors became desirable. 
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&u 
tefls IP that the n functions X, s= ^ form a covariant tensor 

of rank oae; this is known as the towor gradient. If X r is any 
opmiitBt tensor of rank one its simple product by itself or 
" square" is a eovariant tensor of rank two, X n ^ X T X. 



* is an invariant (r, $ umbnJ) 

life Is ctfled the square of the magnitude of the tensor. In 
p*rtJe*rbT the square of the tensor gradient is the invariant 



is known as the " first differential parameter of u." Similarly 
tbe magnitude erf the square of a contravariant tensor of rank 1 
k Hie invariant j Jffr>Z w . 
Again 

^S sZt (rumbral) 

ox 

h contravariant of rank one (rule (d)). Hence 

1 A / fat \ 

^ * ^ an invariant ( r > s umbral ) 



by the result of the preceding paragraph. It is written A 2 w 
and is known as the "second deferential parameter."* In 
ordinary space of three dimensions in which the z's are rec- 
tangular Cartesian coordinates 

g r * = if r 4= s 
= 1 if r = s 

and g n * /*; V^ = 1 so that A 2 w takes the form 



, _ , 

I * /9^9 I 



* Lormor, J., Transactions Cambridge PhiL Soc., Vol. 14, p. 121 (1885), 
obtains tiiis transfonnation in the case n = 3 by the application of the Calculus 
of Variations. 
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When we change over to any "curvilinear" coordinates 
w have under the form 



the expression of this magnitude in a form suited to the new 
coordinates. 

6. GENERAL ORTHOGONAL COORDINATES 

TObenever we have, in any space, coordinates x such that the 
expression (<&)* involves only square terms, Le., g n B if 9 + r, 
the coordinates are said to be orthogonal (for a reason to be 
explained later). It is usual to write, in this case, 

W - A (dxi)2 + Ti (<k (z> )' +--+TI 

fli 2 n? kS 

accordingly 

1 1 



so that 



1 1 1 

t-Ti'Ti-~TT*> 



*. 

= r* 



The square of the gradient is 



whilst the quantity 



_ a |L)} 



The reader should write out the explicit formulae for space polar 
and cylindrical coordinates in ordinary space of three dimensions. 
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(&) F*oaa any sir-vector X n we may derive a four-vector 
alternating tensor of the third rank) 



It is tiiis fosnvvector that was written lor X r * in honor of 



OF THE RECIPROCAL RELATIONSHIP be* 



, 



tween an alternating tensor of rank r and one of opposite kind 
of r*Bk u r from the case r = 1 already treated to a general 
vihieof r. 
We have already seen that 



where J is the determinant 



of the transformation 



J (J, p ) is the cofactor, in the expansion of J, of the element 
(j tf ) == - } of this determinant. 
Hence 

! ^-.CrO *,(n) 

/., J^ 



Now the determinant of the minors of J is well known to be 
equivalent to the product of / by the determinant of order n 2 
obtained by erasing the $ith and $ 2 th rows and the nth and r 2 th 
columns of J affected with its proper sign (the determinant of 



order n 2 is the cofactor of 



in the Laplacian ex- 
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pansion of J in terms of two row determinants from the **th 
and **th rows and the nth and rath rows). Hence we have the 
result that the n(n l)/2 distinct components of an alternating 
covenant tensor of rank n 2 when divided by Vf form the 
distinct components of an alternating contravariant tensor of 

rank two. And so in general. Simikrly the f * j distinct com- 

ponents of an alternating contravariant tensor of rank r 
when multiplied by *>[g form an alternating covariant tensor of 
rankr. 

Example. Take n = 4, r = 2 and consider the linear orthog- 
onal trajisformations of the Special Relativity Theoi-y. Here 



M 



The two tensors or six vectors X r9 and Z r * were said to be 
reciprocal.* 

* Cf . Cunningham, E., The Principle of Relativity, Ch. 8, Camb. Univ. 
Press (1914). 



CHAPTER IV 
L 6iB6S!meAii INTERPRETATION OP THE COMPONENTS OF A 

TENSOR 

DEFINITIONS 

($ Ihiwfam qf a cun% at any point on it. 
Ai any point u on the curve Vi specified by the equations 



whose length 3 from a fixed point 14, is defined by the integral 



s== ' ~!9*n 



du du 
we may form the n quantities 



ZW = = <r . < , , 

* ~~fa"fa (r =1 ) 

We exclude from consideration here the " minimal " curves 
along which d, = 0. Since J* s <faw is a contravariant tensor 
of rank one and ds is an invariant we have that the n quantities 
!, , f rm a rj contr ^ariant tensor of rank one which we call the 
direction tensor of the curve at the point in question. The 

efficients - The 



= ! S that a knowled ^ of ^e mutual 
of the duectam coefficients suffices to determine their 

extV 8 " 6 ^ r indefiniteness to sign). Otherwise 
xpressed the magnitude of the direction tensor is unity. Fixing 
the mdenn.teness as to sign by a p articu l ar choice is J aid tfl fi ^ 
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a ***** f direction on the curve and the curve may be then said 
to be directed. 

2. (6) Metrical Definition of Angle 

Consider two curves with a common point and let their direc- 
tion tensors at this point be P r) and m (r) . The simple product 
Xn ^ P r) i (t) is contravariant of rank two (Rule (e), Ch. 2) rod 
so the expression grJPmW is invariant (r, * nmbral; Bide 
(d), Ch. 2). This we call the-cosine of the angfe 9 between the 
two curves (directed) at the point. If the quadratic differential 
form defining (dsf is supposed to be definite, i.e., if it is supposed 
that (ds) cannot be zero, for real values of the variables x< r) 
and dx (r) save in the trivial case when all the daf^ = 0, it can 
easily be shown that the angle defined in this way is always real 
for real curves. Let us write instead of da: (r) the expression 
XI (r) + ;*m (r) and thus form the quadratic expression in X and /t 



This is not to vanish for real values of X, p save when X = 0, p 
(we suppose the quantities / (r) and m (r) all real and the two direc- 
tions as distinct). Using 

grJW* = 1 = g n mV>mM 
we have that 

X 2 + 2X/i cos B + M 2 = 

must have complex roots when regarded as an equation in 
X : M- Hence 1 cos 2 6 > so that the angle as defined above 
is always real for real directions under the assumption that (ds) 
cannot vanish on a real curve. It must be remembered however 
that this assumption is not always made, e.g., in Relativity 
Theory. 

When cos 6 = the curves are said to be orthogonal or at 
right angles at the point in question. 
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EXAMPLES 

fi* cf&ary space with the ar's as rectangular Cartesian co- 
oadiaates we have the usual expression 



, /*), ( ,*) ^ the direction co^ne, of 
. If BOW we use any " curvilinear " coordinates 
(if f*** J^) the angle between two curves is 



if we me orthogonal coordinates 



Thus for a curve in polar coordinates r, fl, 





It wffl aow be dear why those coordinates in terms of which 
W has no product terms are said to be orthogonal. 
For 



f ma 

9 ^d^faT ^ froin lts co variant character) 



' are 

ve a curve whose equations, in the x coordinates may be 
convemently 8pecified by means Q{ ^ ^ * s > ^ be 



each point , there pass curres of this kind which we 

rf "J?**"!* ^ ^^ that point. On the 

of these coordinate line, the direction tensor 
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aad so the vanishing of the component / states that the co- 
ordinate Unes y< r) and jfi* are orthogonal. Hence if (<&)* does 
not contain any product terms the coordinate lines are everywhere ^ 
all mvtitally orthogonal and so the coordinates are said to be 
orthogonal. In ordinary space, I.e., where the x's are rectangular 
Cartesian coordinates and where the y's are orthogonal co- 
ordinates, 



and 
so that 



a result which is sometimes useful in the calculation of the 
coefficients /n, /22, /ss of the form (<fe) 2 in the curvilinear 
coordinates y. 

3. RESOLUTION OP TENSORS 

If we consider any covariant tensor X r of rank one and take the 
inner product of this into a direction tensor Z (r) we derive the 
invariant Zr/ (r) (r umbral; Rule (d)). This we call the 
resolved part of the covariant tensor along the direction Z (r) . 
Let us now make a transformation of coordinates from x to y 
and consider the coordinate line ^ (a) . The n components of the 
direction tensor for this curve are proportional to 



To determine the actual values of these components we must 
divide through by the positive square root of 



and this is equivalent to 
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equations defining the co variant correspondence for a 
teaeor of Uhe first rank are 

(*= 1> -'^fi; rumbral) 

: ' 

5\ s* $# times the resolved part of the tensor X r along the co- 

r 1 <xrdinate direction y (0 

'- A II 

EXAMPLE 

Space polar coordinates # in ordinary space of three dimensions, 
Tifc s ai rectangular Cartesian coordinates. Denoting the 
resolved parts of the wvariant tensor X in the directions 
f^f J^> y* by fi, 8, $ respectively we have since /n = 1 ; / 22 = r 2 ,- 
/=^sin 1 9 

7i as R; 7 2 = r9; Y 3 = r sin 0*. 

jee distinct components of the alternating co variant tensor 
i two, curl X , in polar coordinates are 



a"r 

Similarly for cylindrical coordinates p, 0, z where /n = 1 ; 
/ = p 2 ; /a* = 1 if we denote the resolved parts of X along the 
three coordinate directions by R, $>, Z we have YI = E; 
7j s p<; 7 8 = Z and the components of the curl are at once 
written down. 

Resolution of Contravariant Tensors. 

To define what is meant by this we require, not as before the 
coordinate lines y (r) along each of which all the coordinates y but 
one, y ir \ are constant, but the coordinate spreads F n _i along each 
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of which all the variables but one, y& say, vary. The parameters 
! * * - _i may here be very conveniently choeea to be the 
coordinates y\ - > #n themselves omitting y, and tiien $* is a 
constant (a particular function of 1*1 tt_i). Now, in general, 
when we have a F_i specified by equations 



fr 1, - -, *> 

we obtain on the spread, through each point, n 1 paiaineter 
lines by letting in turn each parameter vary, keeping all the rest 
fixed. Any one of these, u r varying, say, has at the point in 
question a direction tensor whose components are proportional to 



Let us look for a direction orthogonal at once to the n 1 
directions of these parameter curves. Such a direction tensor 
has components n (1) 7i (n) say and is said to be normal to the 
spread F-i at the point in question. To express the required 
orthogonality we have n 1 equations 

=0 (Z, m umbral; r = 1 - n) 



du T 

homogeneous in the n (1) - n (n) and thus serving to determine 
their mutual ratios. To actually solve divide across by one of 
the unknowns n (n) say and we have n I linear, non-homogene- 
ous equations for the (n 1) unknowns 

n ci> n (r-i) * 



n (n) 

*The algebra following here is somewhat complicated and so it may be 
desirable to derive the expressions for the components of the normal direction 
tensor to the spread y (n) as follows. Working with the coordinates y the n 1 
parameter curves 2/ (8) varying (s = l, -, n 1) have their direction coef- 
ficients proportional to 

(1, 0)1 

(0, 1,0 0)1 

(0, 1, 0) 
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The determinant of the coefficients has as the element in the 
rtfa row aad *th column 



fc-r ( m umbral; r, 9 = 1, , n 1) 

TOr 

determinant is therefore the product of the two matrices 

dx dx> 



fit 



fa 



each f 1 rows and n columns. It is well known that this 
product can also be written as the sum of products of all corre- 

rttpcctiveiy. The - 1 equations expressing that n< r > is orthogonal to these 
1 directioos are 

/<rn (r > =0 (t = 1, - - , n 1; r umbral) 
Hence the ratios 

fltt) ; flCO : - : n^** = /** : /** : : /** 

the actual values being these divided by 



/ n * 
[one most be warned against thinking that ~j= (a 



1, -, n) are contra- 



variant When a change of coordinates from y to s is made the spread 
y<> const does not become x (n > = const.] If now we wish to use x co- 
ordinates, the normal direction tensor, being contravariant of rank one, has 
components proportional to 



; 



ff 



fr - 1, 



umbral) 

If y<*> M y( a; a) | . . . 9 x (n)) w } me fo^ fa 1Mynna i direction tensor to the spread 
ha* it* component* proportional to g* -^ the result required. 
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spending determinants of order n 1 that can be fanned firom 
each matrix. Let us write for brevity 



tia- (-!>*"''*> 



and the determinant of the coefficients becomes 



which may be written f"(JJ . The numerators of the fractions 

furnishing v\ fln-i are dealt with in the same way and we have 

n d) : n <2> ..... n oo B 



(Since the (J,) are really the n distinct components of an alter- 
nating contravariant tensor of rank n 1 we know that 
X* = (J) V^ is a covariant tensor of rank one verifying the 
contravariant character of the n (r) (Rule (d))). If all the (/,) 
vanish the point is said to be a singular point of the spread and 
the determination of n (r) becomes impossible. 

Let us now apply these generalities to the spread V^i given 
by a single equation 

V( x m . . . a-W) = o 

connecting the coordinates x. We may solve for one of the 
coordinates, z (n) say, in terms of the others z (1) - - x (n ~ and 
these others we use as the n 1 independent parameters of the 
spread: 



are then the equations, in parametric form, of the spread PVi. 
Our matrix 

dx (T) 

-(r = 1, ...,n;*= 1, --, n- 1) 
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100 
1 





and 



1 




dx 



But, on differentiating the equation F(a; (1) 
spread Vn^i we obtain 



a:(>) = of our 



dV , dV ds< n) 

\ i \ /..x % M^ 



so that 



whence 



In particular, if the spread FB_I has, in the y coordinates, the 
equation y w = const., we have for its normal direction tensor 



n 







dx (> > 



umbral) 



The actual magnitudes of these components are found by dividing 
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through by the positive square root of 



which expression is = g*' (r TIC* ombral) 



If BOW we have a contravariant tensor Z (r) of rank one it is 
meaningless to call Z (r) Z (r) the resolved part of the tensor In the 
direction I for the simple reason that this expression is not 
iwoariant but takes on different values in different systems of 
coordinates. However, we may first form the covariant tensor 



(r umbra!. Rule (d)) 

This tensor is said to be reciprocal to the contravariant tensor 
with respect to the fundamental metrical quadratic differ- 
ential form and its resolved part in any direction we call the 
resolved part of the contravariant tensor in that direction. Thus, 
for example, the resolved part of the contravariant tensor X T 
in the direction normal to the coordinate spread # (r) = constant is 

93P X n<-> = g^pf*. + Vp (s, p, t umbral) 



} r (r) 



Hence any component 7 (r) of a contravariant tensor of rank 
one is the product by V/ rr of the resolved part of the contravariant 
tensor normal to the coordinate spread i/ (r) = constant. It is 
now apparent that to deal with covariant and contravariant 
tensors of the first rank we require the coordinate lines through 
each point and the normals to the coordinate spreads through 
that point. When the coordinates are orthogonal, and only then, 
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these fiaee and Bonaals coincide and a great simplification is due 
to ifeb fact This explains why orthogonal coordinates have 
used, ablest to the point of excluding all others, in the 
of The(H?etScal Physics. 



(a) 

Space polar coordinates. These being orthogonal the normals 
to the spreads r = const, 9 = const., = constant are the 
coordinate lines r, 0, < respectively and, if we denote the resolved 
parts of He oontravariant tensor Z Ct) in these directions by 
B> 6, f> tbe three components are 



-1 , * . - 

r rsm0 

In general for orthogonal coordinates y with 



we have^f 1 " !//. and if, as usual, we write f rr s l/A r 2 we have 
/=(W AnV 1 and =A r 2 



Here y = ^yy . . . yc) = ^(/y where we denote by 
i) - (fi) the resolved parts of the contravariant tensor along 

the coordinate directions 1, 2, - - - , n respectively. The divergence 

of the contravariant tensor 



takes the form 






Thus, for space polar coordinates, the divergence is 
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sad for cylindrical 



Example (6) 

In order to illustrate the distinction between covariant and 
ortravariant tensors we now consider Mique Cartesian coordi- 
nates y so that 

+ OW + (<V 3) ) 2 + 



where the constants X, p, v are the cosines of the angles between 
the oblique directed axes. Here 



1 v M 
v 1 X 
ft \ 1 



square of volume of unit parallelopiped with 
its edges along the three axes. 



i -c>> ->(/ = sin X cos 61 = sin ju cos 2 ~ sin v cos 3 where 6 1 is the 
angle between the coordinate line y\ and the normal ni to the 
coordinate plane y\ = const, with similar definitions for 2 and 8 Z . 
Hence 




If we have any vector whose components in rectangular Cartesian 
coordinates (z (1) , z (2) , z (3) ) are Zi, Z, J 3 this vector may be 
regarded as either a covariant or contravariant tensor, i.e., 
X l s A r(1) ; X z = Z (2) ; Xz = X& and if we denote the resolved 
parts of this vector along the coordinate lines y by (A^, X^ A\) 
and along the normals of the coordinate planes y by (X ni , A" n ,, 
J n ,) we have 



J T i 
yd) s yjS Z ni s AT ni sec ff x ; F (2) - Jn, sec 2 ; T 3 = ^, sec ^3 
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(F*, Ft, Is) are the reao&wZ parfc of the vector along the 
three coorfimte Haes whilst (F>, F, F (8 >) are tlie components 
tie vector along these same directions. The tensors Y r and 
tie fee^ocul with respect to the differential form (<fo) 2 , Le., 



Y l = F + *F<*> + /iFW, etc. 
Let us now consider the contra variant tensor whose components 

yea H ^yo); F<*> . p 2 y^ ; yw = P3 y(3) 

wbeff |^ pt, A are scalar or invariant numbers; we find for the 
components in the rectangular coordinate system x 



car __ 

Z r = p/Z* (s umbral) 

where 






Now -TTT is a contra variant tensor and -^ is a covariant tensor 
oy w dx ( *> 

if we regard the y's as fixed and consider merely transformations 
on the # J s so that p/, being the sum of three mixed tensors, is 
actually, as the notation implies, a mixed tensor of rank two. 
It was in this geometrical way that Voigt introduced the idea 
which he called a tensor. The mixed tensor pi is completely 
specified by the three directions y and the scalar numbers pi, P2, 
pj. If the mixed tensor is to be symmetric for every choice of 
Pi, PS, PS we must have 



These equations lead to the conclusion that the " axes " y of 
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the tensor are mutually at right angles and so such a tensor was 
called symmetric. 

In order to study the behavior of the vector X as X changes 
direction, keeping its magnitude unaltered, we may solve the 
equations for X and obtain 

X r ^ TT T ~X* 

where from the geometrical construction T/ is a inbred 
with the same axes as p/ but 



, etc., 
Pi 



so that 



Then squaring and adding the equations for X T we find that Z 
traces an ellipsoid, called the first tensor ellipsoid. 

For a symmetric tensor the directions y are orthogonal so 
that YI = r 1 , etc. A simple example of a symmetric tensor is 
furnished by the uniform stretching of a medium along three 
mutually perpendicular directions successively. It was from 
this example that Voigt originally took the name "Tensor." 
Reference may be made to any treatise on the Theory of Elas- 
ticity for an amplification of the remarks of this paragraph. 

5. GENERAL FORM OF GREEN'S FUNDAMENTAL LEMMA 

Starting with any invariant function of position F(x (1) z (n> ) 
we have seen how to form its covariant tensor gradient 

Y dV , 1 

Xr ^d^ (r=l-n) 

the square of whose magnitude is the first differential parameter 
of V 



68 TOBOTOE ANALYSIS AND RELATIVITY 

Now the acraal direction tensor to V(x x (n >) = const lias 
oompooeats whose ratios are 



the actual magnitudes of these being found on division through 
fey tbe positive square root of AiF. Hence the resolved part of 
the cOTmttaat tensor gradient along the normal is 

(MumbJ) 



this is m i(AjF,* This we shall call the normal derivative 
of F and denote by the symbol - The resolved part of the 

gradient along any direction I is 

aT/ 

(r umbral) 



This we denote by -rr and cofi ike directional derivative of F 
ol 

along the direction Z. The angle between n and I is given by 
cos 8 = gJPn<* - -= l^f'ft ^ *> l umbral ) 



= 1 jwlEL = 1 3F 

Hence 

dV dV 4 

s COS u 

showing that the maximum directional derivative is that along 
the normal. (In general, if we say that any covariant tensor X r 
has a direction specified by the reciprocal contravariant tensor 



X 9 = fX T (r umbral) 

* If we define the " direction " of any covariant tensor of rank one as 
that of ita reciprocal contravariant tensor we may say that the gradient of 
any invariant function of position has a direction normal to it. 
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the resolved part of X r along any direction J is the product of tibe 
magnitude of the tensor into the cosine of tibe angle between I 

the direction of the tensor.) 
The contravariant tensor reciprocal to the gradient <rf K is 



Accordingly, on multiplying each of these expressions by ifj, 
we derive the n distinct components of an alternating covenant 
teosoar of rank 1 (cf . Ch. 3, 4) and so we can foam the 



over any spread of n 1 dimensions given by 



the symbol (J r ) denoting as before 
(_ i)~^ 



The normal contravariant tensor to the spread of 1 dimen- 
sions has, as has been shown, components proportional to 

9 r '(J) (r = 1, , n; * umbral) 

the actual magnitudes being found by dividing through by the 
positive square root of 

9im9 l '(J<)9 mt (Jt) (I, , 9, t umbral) 

= 9 8t (J*)(Jt) (9, t umbral) 

/7N dV , dV 

Hence g (/)g^ = product of V^ z (J,)(J)by ^" t" e curec- 



tional derivative V normal to the spread F n _ a over which 7 n _i 
is being extended. Hence we may write 
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vim by rfF^-i we mean the invariant 
&fe~i* <T!it this is invariant follows from rule (d) since V? (J r ) 
fa * oovariant tensor of rank one (cf. Ch. 3, 7).) 
Applying Stokes' Lemma to /n-i we have 



where the integral J-i is extended over any Fn_i which is closed 
aad tbe ict^rai / OE the right is eztended over any region of 
apace F* botaided by F_I. Here 



mad d n is the invariant 

If, instead of the eontravariant tensor X T = g rt ^-^ , we start 
out with 



where Z7 is an invariant function of position we find 



On interchanging the functions U, V and subtracting we have 



which is the usual form of Green's Lemma. The previous 
equation may be written 



V - &(U, V})-dV n 
where A(U, V) is the invariant mixed differential parameter 

j (r, S umbraI) 
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In particular, if the invariant functions U, V are identical we 
have 



last identity is the basis of various uniqueness theorems of 

Theoretical Physics, If we know the values of U or ow a 

en 

doeed F-i as well as the values of &JQ throughout the region 
bounded by Fi*-i the function IT is unique, save possibly to an 
unimportant additive constant. For, applying the last identity 
to the function W ss Ui U\ where Ui and Ut satisfy tibe above 
conditions, we have 



Now under the hypothesis that 



is a definite form we see that &\W is one signed and vanishes only 

dW 

when all r^are zero. Hence since J*&iW-dV n = we must 
dx (T > 

have TT = throughout the region of integration (r = 1 - n), 
dx (r > 

Therefore, IF is a constant and if the values of U are assigned 



on the boundary and so W = or Ui = U^. 

The whole argument depends on the definiteness of (<fe) 2 . 
Suppose we wish to apply the theorem to solutions of the icate 
equation 

^L'j-^!I 7 j- -- == 
dx 2 dy 2 W (f d? 

Here we have 

+ (dy? + (dz? - 
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and so 



and the theorem cannot be applied since AiF can vanish without 
implying the vanishing of all the derivatives. 

6. APPLICATION TO MAXWELL'S EQUATIONS 
One of the most interesting applications of the algebra of 
tensors is the discussion of Maxwell's Electromagnetic Equations. 
These consist of two sets, which in the symbols of restricted 
vector analysis and the units employed by Heaviside are 



(a) - 

c ot 



= 0; divS = 
c ot 

D is the electric displacement, H the magnetic force, and j the 
current vector; B is the magnetic induction, E the electric force 
and p is the volume density of electrification. We take n = 4 
and as coordinates, in the above form, 

X 0) x . 3(2) = y . 3(8) = z . x (4) = t 

If we assume that 

XK = B x ; Xzi = B y ] Xn ^ B z ; Xu = cE x ] X^ s cE y - } 

Xu s cE z 

are the six distinct components of an alternating covariant tensor 
of rank two, the four equations (b) express that 



Z 



124 



-,- 
oy 



-AXti__dB x . fdE v dE x \ 

-r = -TT- + C ( --* -5 J = 0, 

dy dt \ dx dy J 



dt dx dy 

__ &B V . / 

= -- zr+ c (" 
d^ \ 

dE v \ 

-- 5~ I 

dz J 



__ iz .^ -34 _i 41 __ V . j/z * _ n 

134 = -57- + -T --- h ~T = -- zr+ c ("F~~"T~") > 
a^ da: a^ d^ \ dx oz 

dX2&,dX34 , dX&__ dB x . A (dE z dE v 

C 



dz dt \ dy 
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In other words the integral 



, 2) + B*d(z, z) 



is the integral of an exact differential -its value when extended 
over any closed spread V* is identically zero. Hence its value 
vjjen extended over any open spread F" a can be expressed as a 
fine integral J'Xrdx^ round its boundary. Chi writing 



we have 

7 2 /! 

and an application of Stokes' Lemma tells us that 



or 



_-*. O ^y^ar 

"" "' '"" 



The covariant tensor of the first rank (A x , A v , A z , c<f>) is the 
" electromagnetic covariant tensor potential " of which the 
first three components form Maxwell's vector potential, 4> being 
his scalar potential. 

Similarly, if we assume that ( D x , Z) y , D t , cH x , cH yj cH t ) 
are the six distinct functions of an alternating covariant tensor 
XT* of rank two the equations (a) say that 

-3^254 = CJ x ' t Z 3W S CJy] Zi24 s OJ f ] J5Tl23 = P 
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and we have J 2 = I* where 

It - ya*fc t) + cHJ(y, f) + cHJ,(z, t) - 



, x) - 
, y, f) p d( 



Is bring taken over any closed spread F 2 of two dimensions and 
I, bring taken over the open F 8 bounded by F 2 . According^ 
(^ ^ jV, - P/C) are the four distinct functions of an alternating 
tensor of rank three and so, on writing c(Xi) === jr 23 4, etc, ^ r l 



form a contravariant tensor of rank one (Ch. 3, 7). F rom it$ 
^Mnition we know that its divergence is zero. This tensor we 
may call the current contravariant tensor and write 



.uet us now apply these methods to the problem of writing 
Maxwell's equations in a form suitable for work with curvilinear 
coordinates y a) , y\ y in space of three dimensions the time t 
not entering into the transformation. The equations connecting 
the x and y coordinates are of the type 

Z< 2 >== xM(yW 9 y&> 9 y<3>). 

and denoting tensor components in the new coordinate system 
by primes we have 

(Bo'-ft,'-ft,?gr 



the terms in /T 1} F 2 , T 3 vanishing since 
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the terms in (Pi) (D) (#3) vanishing since 



Hence in the three-dimensional space with coordinate systems 
(a, s, cs) ) and (y (1) , 2/ (2 >, y) the variable I bring regarded 
merely as a parameter which does not enter into 



* fr> * = 1, 2, 3) 

the three quantities (Di) (1)2) (#3) are the three distinct members 
of an alternating covariant tensor of rank two. Hence ~ s= X T 

/ jj \ _ 

is a contravariant tensor of rank one; similarly ^4- = X r is a 

V^ _ 

contravariant tensor of rank one whilst E r = X T and H T = X T 
are covariant tensors of rank one. We derive by our rule (d) 
of composition the invariants 

(ED) (EB). 

V^ J V^ ' V 

where as in the usual vector notation 
(ED) s E^ + 2 Z) 

and similarly for the others. 

Dividing Maxwell's equations, as usually written, across by 
V^ we obtain 

- 1 * Y r + 1 curir (S) = C' (r = 1, 2, 3) 
cdt ijg 

(where C r ss ^is the contravariant current vector). 
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div X (r > p 
p b the invariant charge density and similarly from the 



second set 



_ 
div X r = 

Denoting, then, as usual resolved parts along the coordinate lines 
by subscripts (/i, fe, W and along the normals to the coordinate 
surfaces by the subscripts (n\, n*, w 3 ) we have the three equations 



The equation div Z r p becomes 



(by 1)^ is meant the resolved part of the contravariant tensor 
D/V<? along the direction ni). 

The equations (t) are similar and are simplified by the fact 
that there C ni , Cn,, C ni , p are all zero.* 

* When the coordinates y are orthogonal 



; f 1 = Ai*, etc., and Maxwell's equations become since ni *= li, etc, 
d Hl 



and two fdmilar equations together with 

d f Dl i\_L. d / D| \_i_ 3 / D <3\1 

^ v/ai; + -jja (jtfj + en (ffc) } -p 



Q. -RieTnann-TTefrer, Die Partietten Dzfferentialgl&ichungen der Mathemat- 
itchen Physik, Bd. 2, p. 312 (Vieweg & Sohn) (1919). 
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EXAMPLE 
In space polar coordinates Maxwell's equations are 




It is particularly to be noticed that Maxwell's Equations are 
essentially of a non-metrical character. No real simplification 
is introduced by the hypothesis that the fundamental space is 
of the ordinary Euclidean character. Another point to which 
attention should be directed is the difference in character of the 
tensors B and H or of D and E. A relation of the familiar type 



p, the coefficient of permeability, being supposed invariant is 
not the proper mode of statement of a physical law if we under- 
stand by Bi, J5 2 , J5 3 the three components of the tensor B. The 
true statement of the law is 



where by (B)i we mean the resolved part of the contra variant 
tensor (B)/V0 along the direction I and by (ff)i we mean the 
resolved part of the co variant tensor H along the same direction. 
Thus any constitutive equation of this type is an allowable state- 
ment of a physical law not because it is a tensor equation (since 
it is not such), but because it is an equality between invariant 
magnitudes or a scalar equation. The true tensor equation is 
found by equating the covariant tensor pH to the covariant 
tensor reciprocal to the contra variant tensor 



CHAPTER V 
L COHHCTO>K OF TBKBOB ALGEBEA WITH INTEGBAL 

EIAJNTS AK3> APPLICATION TO THE STATEMENT OF 

FAEABAY^S LAW OP MOVING CIBCUITS* 
Suppose for cacainple we have a curve V\ whose equations 

r) (5= 1, -..,n) 



& porama&r r. This curve may be said, adopting the 
language of dynamics, to move and trace out a V z whose equa- 
tions are those given above, the parameters being u and r. 
Any one of the curves r = constant will then be a position of the 
moving curve. We shall suppose that the values of u serve to 
identify the various points on the moving curve; thus if u 
denotes the distance along the initial position of the moving 
curve from a certain fixed point, or origin, the curves V\ obtained 
by taking u = constant (UQ) in the equations 

*<*> - z ( *>(>, r) (* = 1, ..-,n) 

are the path curves of the definite point on the curve I 7 i which 
initially was at the distance UQ from the origin on Vi. It will 
fix our ideas to Consider Vi as made up of particles of a fluid; 
then the curves V\ are the paths of the various material particles 
of Vi. It is well to insist, at the outset, on the point that the 
parameters u and ? are independent. Thus if the moving curve 
l\ were rigid, u could be taken as the arc distance along V\ at 

* An elementary presentation of the theory of Integral Invariants is given 
by Goursat, E., in two papers: 

(a) Sur les invariants int6graux. Journal de Mathematiques, 6 e sne, 
t. IV (1908), p. 331. 

(6) Sur quelques points de la th6orie des invariants intSgraux. Journal de 
mathe'matiques, 7 s^rie, t. 1 (1915), p. 241. 
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any time r ; if, however, as in tibe case of the curra made op of 
uiaterial jflmd particles, FI is not rigid, u may only be taken as 
tlie initial identifying are distance; otherwise tt would way with 
r. Let us now consider an integral Ii ss fX f dsf^ extended 
over Vi and ask the conditions that Ji should be the same for 
all the curves FI, i.e., that I\ should not vary with r. If this is 
so, h is said to be an integral invariard. 
Now Ii is in general a function of r defined by 



Ii(r) s J 



du (r umhrai) 



the limits UQ and u' being, however, since u and r are independent, 
quite independent of r. Hence 



dr UQ dr\ du 

The coefficients X T are functions of position and therefore involve 
r indirectly; it is somewhat more general to contemplate the 
possibility that they may involve r, not only in this indirect 
manner but also directly. Then for any one of the coefficients 
X r we have 



7 *" 1 *v I *\ f \ *\ 

ar \or dx w or 

It is now convenient to denote the contravariant tensor of rank 

J /y,\ 

one by the symbol X r and to use the result 

dr 



dd^ 8 dx^ ( , n 

-- = -- = - = -T -- (s umbral) 
dr du du dr du dx ( <> du 

and we have 



dr \ dr du dr du 



r 

du r dx^ du 

(r, s umbral) 
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(on modifying suitably the umbral symbols) 

Heoee if dJi/<ir ia to vanisli identically for all curves Five must 
lams 



t it is only necessary that Ji should be unchanged for 
aB ofo^rf curves Fi; in this case I\ is said to be a relative integral 
invariant. To find the conditions for this we use Stokes 9 Lemma 
to replace the Ji over a dosed curve by an J 2 over an open Ft 
and then find the conditions that J 2 should be an (absolute) 
integral invariant. 
The analysis necessary to find the conditions that an 



extended over a F p (moving) whose equations are 

*<> K xM( Ul - - tip, T) (* - 1 n) 

shoidd be an absolute invariant is identical with that given for 
the simplest case p = 1. Let us write as before 



_ = 

ar 

and denote by the symbol F the derivative 
dF_dF, dF 



where F is any function of position which may also involve r 
explicitly. Then 

.. umbral) 
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gfoce the limits of integration with respect to the variables t* are 
independent of T. This we write 



d 



.,..^ 

and availing ourselves of the relation 
d dxM 9 



dr 
we arrive at the conditions expressed in the form that 

f 

A--* -1- A^- 



- = (m umbral) 



An especially simple case is that in which p = n. Here there is 
a single condition 



SdXV\ , . lx 

i... n f ^- J (r umbral) 



Since JTi... n is the single distinct member of an alternating co- 
variant tensor of rank n 



where U is an invariant function of position and writing out 

i... n , 
-- r 



i... n - - * ,,\ 
dr dx (t > 

our condition that J*U-dV n should be an integral invariant may 
be written in the form 



or on dividing out by V<7, which does not involve r explicitly, 
r h div (C7Jf (r) ) where as usual the divergence of the contra- 

C7T 
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variant tensor of rankoiae UT r is the invariant -=^ (V UX*). 



Jm this form the invarianee of the condition for an integral 
invariant & apparent. If we are considering a moving charged 
mteriai body where p is the density of charge, the total charge 
remaining conrfatii gives us that 



at 
where JT (r) is tfae contravariant velocity tensor of rank one. 

Far&day'* Law for a Moving Circuit. 
We have seen that 



the integral in each case being taken over the position of the 
moving curve at time T. The expressions 

BX r , t dX r , 



must accordingly form a covariant tensor of rank one. In fact 
we may write this as 

d3>r , Vt \dX r dX 8 I . d , v v(8 )) 

= a7 +A 9^-) + { ^ Y()! 



when the covariant character is apparent by rule (d), Ch. 2, 4, 
since 



is covariant of rank two and A". A" (a) is invariant. 

Let us now write down the expression for dl^/dr where J 2 is 
any surface integral and transform the coefficients as above so 
as to make evident their tensor character. 
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Writing 
eget 



where 



tfee integrals being in each case extended over the position of tte 
moving spread or surface V* at time r. We may write 



yco -r* 4. *t i r 

^^ 



where we have availed ourselves of the alternating character of 
X r t- The covariant character of X T9 then follows from rule (rf). 
We shall apply this result to the surface integral 



n = 3 

so that (Z>i), (2)2) > (^3) are the three distinct members of an 
alternating covariant tensor of rank two. Hence D (r} = (D r ) V </ 
is a contra variant tensor of rank one. The co variant tensor of 
rank one whose curl appears in the expression for X TS is 

A" r mA rm (m umbral) 

so that its first component is 

A r i*Y (2) + A r 13 A r(3) = -Y ( (D a ) - X^(D*) 



It accordingly appears as that derived from the outer product 
of the velocity contra variant tensor and the displacement contra- 
variant tensor. 
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Tbe expression 



If now we o***iw as Maxwell's equations for the moving material 
caedium 

i(Z)) = 
at 

where ( j) is the alternating covaxiant current tensor of rank two, 
so that ( j)/Vis the eontravariant current tensor of rank one C*, 
we have for % n the equations 



Vf + Z V? P + { 



etc. 

Using Stokes' Lemma to transform the surface integral of the 
part in face brackets into a line integral as well as that involving 
curl H in dD^/di we find 



, a:')] 

The integrand in the surface integral on the right is found by 
writing r, 5, t in cyclic order and summing the terms corresponding 
to r = 1, 2, 3 respectively. (The line integral is to be taken over 
the boundary of the moving surface.) The eontravariant tensor 
pA r(r) is called the convection current. In exactly the same way 
we obtain, on making a similar assumption as to what Maxwells 
equations should be for moving media, 



+ 
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there being now, however, no surface integral on the right-hand 
side. Accordingly the covariant tensor 

E r + - V0(-^ (t) B (l) ~" -X^B^) (r = 1, 2, 3; r, *, t cyclic) 

is taken as the effective electric intensity along the Bjoving curve; 
its line integral being called the effective electromotive force 
round the curve. (X {r ^ is the contravariant velocity tensor.) 
On multiplication by charge this tensor gives the mechanical force 
tensor. 

Example. In space polar coordinates the mechanical force 
tensor per element of length on a moving curve with linear density 
<ris 



<rds 



I rE e + -r 2 sin e(v+B r - v r BJ -4-r \ 
[ c r sin0 ] 

I r sin OE+ + - r 2 sin 0(v r B e - v 9 B r ) - \ vds 

where E Ti B r , ^ T are the resolved parts of E, B, X along the direc- 
tion r and so on. To get the resolved parts of the mechanical 
force along the three coordinate directions multiply these by 1, 

-, : respectively and we obtain the well-known formula 
r r sm 6 



= E+-[vB] 
c 



In the general case when the coordinates 2/1, 2/2, yz are not orthog- 
onal the three resolved parts of the mechanical intensity (covar- 
iant) tensor along the coordinate lines y\, y^ yz respectively are 



F h ^E h + - - K n , - fln,BnJ> etc., 

c * ju 

where v ni v^ z? n , denote the resolved parts of the velocity along 
the normals to the coordinate surfaces y\ = const., y z = const., 
2/3 = const., respectively. 



CHAPTER VI 
L THB TENBOB OB ABSOLUTE DIFFERENTIAL CALCULUS 

Since the Calculus of Variations deals with properties of curves 
and surfaces without making any particular reference to the 
special coordinates used in describing the curves there must be 
a doee relationship between that subject and that which we 
ftfe discussing. It is this absolute or tensor character of the 
calculus of variations that has urged writers on Theoretical 
Physics to express the laws of physics, as far as possible, in the 
language of the Calculus of Variations. However, this subject 
has been placed on a clear and firm basis only within the past few 
decades and so it may be well to discuss one of its simpler prob- 
lems the more so as the solution of this problem is involved in 
the statement of Einstein's fundamental law of Inertia in the 
Theory of Relativity. 

Let us consider a curve Vi, in space S n of n dimensions, given 
by the equations 

( S = 1 ... n ) 



and in connection with this curve a function, not merely of 
position, but of the coordinates x and their derivatives 

<*' dx(9} 

() =E _ 

du 
The integral I\ over the curve V\ where 



has a value depending on the curve J"i as well as on the particular 
function. The problem we wish to discuss is: What, if any, are 
the curves V\ making, for a given function F, I\ a minimum, all 
the curves Vi being supposed to have the same end points. 

86 



THE ABSOLUTE DIFFERENTIAL CAJLCDLD8 8? 

To answer this question we consider a aew carro *(a) gran Iqr 

the equations 

a) ( S =I... H ) 



where a is quite independent of u. We suppose 
a to be such that when a = 0, FI(<*) makes /i a minimum. 
Vi(a) is now completely determined by the equations just wtittm 
when a is given and so Ii is a function of (a) which may, we cop- 
pose, be expanded by Taylor's Theorem in the form 



This is written 



and 81 is called the first variation of the integral. If 1(0) is to be 
a minimum it is necessary (although not always sufficient) that 
57 = for otherwise A7 ss I(a) 7(0) would change sign with 
a when a. is sufficiently small. Now the limits of the integral 
for /i are fixed and so to find d I [da we have merely to differen- 
tiate the integrand F with respect to ex. F involves a, not 
directly, but indirectly through the coordinates x and their 
derivatives x r . 
Thus 

dF dF dx^ , dF dx<*' , , n 

T-= T-77)^~ + ~TT7~5 (jumbral) 

da dx (8 > da d x (*r da 

and therefore 

a/ = r u 

^0: J uo 
Now 



da dadu du da 
so that, on integration by parts, 



r F ^ 

J.i dx (>)> da 
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SiBoe the end poiats of the curve are fixed, dx^/da = at the 
Emits of integration and so the integrated part vanishes and, 
collecting terms, we have 



A is to be aero for oH possible varied curves F(a) it is 
sufficient and can be shown to be necessary that all the 



Ct-ip >)< 

These expressions are the components of a covariant tensor 
of rank one where now, however, the term is used in a wider 
sense than hitherto, F is now not merely a function of the 
coordinates x but of their derivatives x'. Prom 



we have 

dx"' 



so that 



dF 



n ri 

showing that : = X r is a covariant tensor of rank one. 

dx (r) 

Suppose we wish to find the geodesies of our metrical space S n . 
These are the curves for which the first variation of the length 
integral is zero. 



the gi being functions of position. We shall find it convenient 
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to take as parameter u the arc distance * aloag the doogfat-lor 
geodesic.* Then when we put a = after the 
F = 1, from the definition of arc distance s, and we 



where we use x to denote ( -r- J so that Z* a (i) tfae 

V < /^> 

unit contravariant direction tensor along the soogbt-foc geodesic. 
Also 



o 
and our equations are 



or 



(r, m umbral; i = 1 ) 

Multiply through by g pt and use t as an umbral symbol so as to 
obtain the n components of a contravariant tensor of rank one 

Wi ^_l|g=) = (r, m ,t umbral) 



It is now convenient to introduce the Christoffel three-index 
symbols of the first and second kinds defined as follows: 



(b) {rs, t] = {, <} = g tp [rs, p] 

* However, this rules out those minimal geodesies along which a is constant 
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nMefa equations imply 

**{**, *} s gtdf*ir*, P] ft P umbra!) 

- fa 0) 
Equations (a) give 



Iben we may write 



[m, t] - |[if, m] - |[te, r]] 
m] + [rm, fl - |M, r]] 



since an interchange of the umbral symbols r t m in the last 
three-index symbol leaves the summation unaltered. 

Accordingly, on using the definition (6), the differential equa- 
tions of the geodesies are 

Z ( P> + {rm, p}dfx m = (p = 1 n) 

From their derivation we know that these equations are contra- 
variant of rank one. We proceed now to obtain a general rule 
which makes the tensor character of equations of this type 
apparent on inspection. 

2, THE FORMULA FOR COVARIANT DIFFERENTIATION 
From the covariant character of the g rs we have 

(I, m umbral) 

, ., *U) ;u(m) 
+ gim ( _^J 



"() 5~()J~(T) ) ^ m> n um ^ ra ') 
where in the differentiation we have remembered that gi m is a 
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function of the #'s only indirectly through the ar'a. We easily 
obtain two other similar equations by merely interchanging 
(r, t) and (s, t) in turn. We are careful to so distribute the 
umbral symbols Z, m> n as to facilitate combination of the three 
equations obtained in this way. Thus 




Now adding the first two of the equations and subtracting the 
third we have, on writing 




Z, m, n uinbral] 
Now 



from its covariant character (p, q umbral) 



To remove the coefficient of _ . . n ... multiply across bv/ jt 

dy (r) dy (t > 

and make s and t umbral when we get 



from the relation (contravariant) 
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Knafly 



from which on interchanging the r61e of the x and y coordinates 
we have 



Seppoee now we have a covariant tensor of rank one X T so that 

v - 

' = 

Hiea 



dY ' 



on altering suitably the umbral symbols Im to rp. These 
equations state that 

-- X k {rp,k} s X rp 



is a covariant tensor of rank two. Consider now a contravariant 
tensor of rank one so that 

^* s ^ f d~lr) ( r um bral) 

Then 



(r, p umbral) 
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These equations state that 



is a mixed tensor of rank two. 

These tensors of rank two are called the 
of the covariant and contravariant tensors X r and X* respectredy. 
Similar analysis can be carried out to obtain the cowiattt 
derivative of a tensor of any rank and character. To make this 
perfectly clear let us take the case of a mixed tensor JT/ of rank 
two: 



. 



whence 



_ , F p( , , , 
= ' ~ Xk{ql ' k] 



expressing that ^fr- - X k v {ql, k} + X q k {kl, p} is a mixed 
ox 

tensor of rank three being covariant of rank two and contra- 
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variant of rank one. In general, the covariant derivative of 



It wiH be noticed that + signs go with the contravariant symbols 
and negative with the covariant. Also the new label s is always 
seoood in the three-index symbols; the umbral label is first if 
taken from the oontravariant and third if taken from the co- 
vari&nt indices. 



3, Amjomofls OF THE EXILE OF COVARIANT DIFFERENTIATION 

(a) Ri&nann's four-index symlbok and Einstein's Gravitational 

Tensor 

Prom any covariant tensor X r we obtain as its covariant 
derivative 

a xr 

X n = gl - X k {rs, Je\ (k umbral) 

and as its second covariant derivative 



t> Pl ~ Xr * {st ' P} 



r) V 

~> 



From this by the elementary rule (6), CL 2, 2, of tensor 
algebra we derive a new covariant tensor Y ra/ = X rt9 and the 
difference of these is a covariant tensor of rank 3 by rule (a), 
Ch.2,l; i.e., 

"Y = i r Y 

A>rst A T9 t -Ar< 
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the terms involving the derivatives of the X T cancelling 
pfetely out. Now X* is an arbitrary covariant tensor 0f rtak 
one and so by the rule (e\ Ch. 2, 5 the converse of the nde 
(d) of composition 



is a mixed tensor of rank four of the type indicated by the 
positions of the labels. 

If we write k = t and use t as an umbral symbol we derive by 
rule (d) Einstein's gravitational covariant tensor of rank two 

G r9 E= A. {rt, t} - Aj { rs , t} + {ps, t} [rt, p} - \pt, t\ \n, p\ 

The mixed tensor X T8t k is usually denoted by the symbol {rt, is] 
and is known as the Riemann four-index symbol of the second 
kind. From it we obtain by the rule of composition the co- 
variant tensor of rank four 

[rj, ts] = g 3 kXr*t k = g,k{rk t is] (k umbral) 

which is known as the Riemann four-index symbol of the first 
kind. From Einstein's tensor of rank two we obtain the in- 

variant 

G == g rs G r8 (r, s umbral) 

which has been called the Gaussian or total curvature of the space. 
This name is given since G is regarded as a generalization of the 
expression given by Gauss for the curvature of a surface ij.e., 
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n 2), The tenn curvature is widely used in the literature of 
aod so it may be well, in order to avoid a possible 
of ideas on the subject, to discuss briefly what is 



meant by the curvature of a metrical space. To do this it is 
necessary to say a few words about the four-index symbols. 
We haw, by ctefinition, 

' k} 



+ {pr,t}{ts,k}- 
Recalling that 
we have 

& (iumbral) 



from definition of [g*, k] so that on operating similarly with 
g^-2-- [ps, k] and writing g qk {ts, k] = [is, q] we find 



(the terms {pr, t\[ts, q] and {pr, k}[k$, q] cancel since i and k 
are merely umbral symbols) . Finally, in terms of the three-index 
symbols of the first kind, 
3 ^ 

lW> 1 = ^(7, t? r . ?! - few ^ 5 ' * ] 

+ ? tj ([F, jlfor, i] - [pr, j][?, i]) (fc, ; umbral) 
Writing out, in the first two terms of this expression, the values 
of the symbols, e.g., 

r i I d d d 

lpr> 9l = 
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we find 



From this formula it is apparent that 

(a) An interchange of the indices or labels p, q 
tie sign of the symbol. 

\pq, r$] + [qp, rs] s o 

(fc) Similarly 

\P9> **] +\P9>*r] = 

(c) A complete reversal of the order of the labels does not 
alter the symbol \pq, rs] = [w, pg]. This depends on the sym- 
metry relations g kj ' = g jk . 

(d) If we keep the first label fixed and permute the other three 
<%dicatty we get 3 symbols whose sum is identically zero, i. e., 

[p?, rs] + [pr, sq] + [ps, qr] s 

The number of non-vanishing symbols which are linearly distinct 
now follows. If p = q or r = s the symbol vanishes on account 
of (a) and (6). The number of choices for the first pair (p, q) is 

nz = ^r - and similarly for the second pair (r, s). However 

L 

relation (c) shows us that we do not get n 2 2 symbols by combining 
the two choices but 



The relation (d) will still further reduce the number of linearly 
distinct symbols. When the indices or labels p, q, r, s have 
numerical values which are not all distinct the relation (d) merely 
reduces to a combination of the relations (a), (6), (c). There are 
therefore n(n l)(n 2)(n 3) new relations in (d). How- 
ever since there are three letters q, r, s permuted cyclically, each 
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relation will occur three times. Each of the relations (a), (6), (c) 

reduces tibe number r which remains in 

o 

half and so there are 

m m ( !)( 2)(n 3) -5- 24 distinct relations (<i). 
; accordingly but 



distinct Riemann fouwndex symbols. For n 2 there is out 
one which we may write [12, 12]. When we change the coordi- 
nates firom x to y we have 

no ,01 r 
[12, 12], - fa, 



(from covariant character) 

Snce there is but one distinct symbol [pq, rs] it will factor out 
on the right and we get (since there are but four of the symbols 
which do not vanish) 

[12, 12], = [12, 12]- J 2 



where J is the Jacobian n , ; We have already seen that 



ri2 121 

/ = gJ* and on division we obtain the invariant K = = - - 

9 

It is this invariant which Gauss called the total curvature of the 
space of two dimensions under discussion. 
In order to compare this with the invariant 

9 r *G n (r, $ umbral; n = 2) 

we have 

G ra = {rt,t8\ (t umbral) 



.: G'n = ^[1221] 

{since if p = 1 or t = 1, [Iptl] = by relations (a) and (&),} 
= - <7n[12, 12] - , 
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feosa definition of p 22 , 

= - gn-K 
Similarly 



12 = <?*[&> 12] = - frrK = - g 
<?2i = 0[21, 21] = - ja-Z = - 

from relation (c), 

<?22 - ^[21, 12] = - <te-K 

so that 

" = - Kf* Sn = - 2K 

since 



= - 2K (r ombcal) 

For a space in which, in some particular coordinate system , 
the coefficients g r * are constants all the three-index symbob 
[pr, $] and in consequence all the symbols {pr, 9} and also the 
four-index symbols \pq, rs] and {pq, rs] = 0. On account of the 
tensor character of these latter symbols we know that the Rie- 
mann tensors \pq, rs] {pq, rs} will be zero no matter what the 
coordinates are. Conversely the vanishing of the tensor [pq, rs] 
expresses the fact that it is possible to find coordinates y such 
that the f r a defined by the equations 



shall be constants. We may now apply the well-known method 
of reduction of a quadratic expression to a sum of squares (as in 
the determination of normal vibrations in dynamics where the 
expression for the kinetic energy is reduced to a sum of square 
terms) ; the transformations on the r/'s are linear in this operation 
and we finally get 



(If we restrict ourselves to real transformations there may be some 
negative squares; thus in the relativity theory there are three - 
and one + term.) A space of this character is said to be Ewlidean 
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tbe #'s are called orthogonal Cartesian coordinates. Rie- 

B defines curvature by means of his tensor [pq, rs]. When 

tbfa tensor radishes the curvature of the space is said to be zero 
do that Eudidean space if one of Zero Riemann Curvature and 
oomertdy. If tl*e ratio of the component [pq } rs] of the curvature 

tensor to the two-rowed determinant 9pr gp * is the same for 

aB values of p, g, r, *, Eiemann says the space is of constant 
otherwise the curvature will be different for different 



orientations at a point:* Gauss' total curvature, on the other 
hand, has a numerical value at each point in space and has 
to do with the different orientations at that point. We 
som up by saying that a gravitational space is, at points 
free from matter, non-Euclidean, i.e., it has a Riemann curvature 
but its Gaussian curvature is zero. 

It may be well to call attention to the fact that the definition 
* The differential equations of the non-minimal geodesies of any space are 



* being the arc length along the geodesic. It is known that the solutions x< r > 
of these equations are completely determined by the values of x (r > and -^~ 

for a particular value of s, * - let us say. This is stated geometrically by 
saying that through any point in space there passes a unique geodesic with a 
given direction. If, now, through a definite point we construct the geodesies 
with the distinct directions (r) and 7j (r) respectively (r = 1, , n) and con- 
sider the family of geodesies through the point in question obtained by assign- 
ing to each a direction tensor whose rth component is proportional to 
X (r) -f M^ (r> and then letting the ratio X : /t vary, we obtain a geodesic spread 
T r j of two dimensions which at the point in question has the orientation deter- 
mined by the two directions and 17 through the point. It is the curvature of 
this geodesic Vi that Riemann calls the curvature of the space relative to the 
orientation determined by and 77. There is a remarkable theorem due to 
Schur i Math. Ammlen, Bd. 27, p. 563, 1880) which says that if at every 
point the Riemann curvature of space is independent of the orientation the 
curvature at all points is the same. Such a space is, then, properlv called 
a spnre of constant cuivaturc. 



THE ABSOLUTE DIFFEEENTUL CALCUI4B 101 



space given above is a " differential " deiraitioci; 
aces which are Euclidean according to this definition do not 
jcessarily satisfy the postulate that one can proceed indefinitely 
a given direction without coming back to the starting point, 
be simplest example is the well-known one of a cylinder dF unit 
dius. In this case n = 2, y a) = <j>, the longitudinal angle 
easuied in radians, and y s 2, the distance measured parallel 
> the axes of the cylinder: 



1 the cylinder is cut along a generator and developed on a plane 
will cover a strip of breadth 2v on the plane. If we take 
octangular Cartesian axes in the plane, with the zf axis parallel 
> the strip, points whose z (2) = <t> differ by 2r correspond to & 
nique point in the strip (that one with the same a;^) and to a 
nique point on the cylinder. Hence there are an infinity of 
traight lines (i.e., geodesies) joining any two points (with 
iff erent z's) on the cylinder. They develop into the oo l straight 
nes joining the points 

(x, z< 2 >) and (*<, z (2) + 2ri7r) ( = 1, 2, - -) 

n the plane. It is evident that speculations as to the " finite- 
LCSS " or " infiniteness )J of a space based on its di/erential 
haracteristics must be regarded with distrust. 



CHAPTER VII 

1. In this final chapter we shall treat in a brief way, as an 
application of the preceding analysis, the classical problems of 
Relativity. As in other applications of the methods of mathe- 
matical analysis to problems in physics the first, and here the 
most serious, difficulty is that of giving a physical significance 
to the coordinates. All systems of coordinates are, without 
doubt, equally valid for the statement of the laws of physics 
but not all are equally convenient. It is reasonable to suppose 
that for a given observer of phenomena a certain coordinate 
system may have a direct and simple relationship to the measure- 
ments he makes; such a coordinate system is called a natural 
system for that observer. It is necessary to define this natural 
system and to find by experience, or otherwise,* how the natural 
systems of different observers are related. This has been well 
done in the special or " Restricted Relativity Theory " but in 
the more general theory, which we propose to discuss here, 
much remains to be done in this part of the subject. In what 
follows we shall consider (a) the problem of determining the 
metrical character of the space-time continuum round a single 
gravitating center and (6) in consequence of the results of (a) 
the nature of the paths of a material particle and of a light ray 
in a gravitational field. We shall, following Einstein, make the 
fundamental assumption that the space which has a physical 
meaning or reality, i.e., with reference to which the laws of 
physics must have the tensor form (cf. Ch. 2, 1), is one of four 
dimensions (commonly referred to as the Space-Time continuum). 

*The relationship between the different systems may be arrived at by 
making various hypotheses whose truth or falsity must then be tested in the 
light of experience. 
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2. THE METRICAL SPACE ATTACHED TO A SBTCOJR GRAVITATING 

CENTEB 

We assume that for an observer attached to the gravitating 
center one of the four coordinates, af& say, of Ms natural system 
is such that the coefficients gu, g&, g& of the quadratic differential 
form for (&) 2 vanish identically whilst those remaining ore 
independent of z (4) ; # (4) is said to be a time coordinate aad the 
field is said to be statical. Accordingly 



ft m - 1, 2, 3 

Now in any space of three dimensions we can always find orthog- 
onal coordinate systems; for, writing the metrical (d*f in its 
reciprocal form (<fe) 2 = P^rty, we have merely three equations 
J** = (r =t= *) or explicitly 



* -^- U 



^to determine the three unknown functions y of x so that the 
coordinate curves y may be orthogonal. There is no lack of 
generality, then, in writing (&) 2 for the statical field in the 
orthogonal form 



where we have dropped the double labeling as unnecessary. (In 
general it is impossible to find orthogonal coordinates in space of 
four dimensions since there are now six differential equations 
fr* =0 (r + s) for the four unknown functions y and these 
equations are not always consistent.) We must now go through 
the details of evaluating Einstein's gravitational tensor (cf. Ch. 
6, 3) for an orthogonal space. 

The relations g r * = 0; g r * = if r =|= s make matters com- 
paratively simple. We shall use r, s, t to denote distinct numer- 
ical values of the labels. Then 
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[rt, t] & ^[rs, Jc\ m g"[rs, t] = Q (k umbra!) 
{, r \ m {-, r\ m j*fr, ft] - fl*. '1 = 



(i being the only umbral label here) 
i t _ J d *' 

{rr> r} ~ 255 

simikriy and 



Riemann four-index symbol of the second kind (cf . Ch. 6, 
3) is defined by 



to* n} m a> {pr > q} ~& {ps ' ql + {pr * l] {h> q] 

{ps, 1} {lr, q} (I umbral) 

and those components vanish identically, for an orthogonal 
coordinate system, where the pq, rs are distinct; [ {pr, 1} vanishes 
unless I = p or r in which case {Is, q} vanishes]. To evaluate 
the remaining symbols write r = q without, for the present, 
using q as an umbral symbol 



The formulae from this on take a simpler form if we use the 
symbols H defined by g r = H r 2 ; thus 

1 PH q __ 1 6H P dH Q 

1 dH 8 dH q 
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we find 

, , H p &H f _ H P dH f dS, 



_ j^ &H, , H p &H f _ H P 
{pq ' W ~ Hi'dx^* H*dxt H, 

q , E, \ 1 9H, 9H t 
"*" H, [H f *dxPds< 



dH f 



where r and are the two labels different from p aad q. Tbe 
components of the Einstein tensor are now found by gmmnmg 
with respect to q. It will be recalled that fep, r*\ (p, r, * 
any values distinct or not, cf. Ch. 6, 3). Hence 



, rs\ m g* k [pk, n] s ywfep, r] = 

SimUarly {^>g, 55} = 0, so that in forming <?, for example, 
we need merely write 

<? 12 = {13,32}+ {14,42} 

whilst 

G u = {12,21} + {13,31}+ {14,41} 

It will be observed that differentiation with respect to x (?) and 
z<> occurs in every term of {pq, qs] and so the absence of the time 
coordinate z (4) from the coefficients makes GH, G&, G all iden- 

tically zero. 

We shall now make the following hypotheses of symmetry (a) 
we shall suppose that the coordinate lines * (1) are geodesies of the 
space (all passing through the gravitating center). The equa- 
tions of the non-singular geodesies have been found to be (Ch. 6, 

D 

(D + {i m> rlx^x^ = (r = 1, - , 4; /, m umbral) 

where dots denote differentiations with respect to the arc distance 
which we take as our coordinate z (1) . Writing 
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(since $f$, a^ f ^ are constant along the coordinate lines 
we find { 11, r J = which from the values given for this symbol 
yields #L = constant Hie constant is in fact unity since, 
by hypothesis, di = dar a) along the curves x = const., 
$* * const., ^ const. It is apparent that it is sufficient 
that 0i be a function of a^ alone for we may make a change of 
variable d& ss zftofy) leaving the other coordinates unaltered; 
the argument shows conversely that if ^i is a function of x^ 
alone the coordinate lines z (1) are geodesies, the arc length along 
them bring given by 9 = /V0T dx. 

{&) 3/8 and a^ are directional coordinates serving to locate a 
point on tbe geodesic surface aft = const., a: <4) = constant. We 
shall suppose that the arc differential on this surface (which 
may conveniently be called a geodesic sphere) cannot involve the 
" longitude " coordinate x (i) nor can the arc differential along a 
given " meridian " z (3) = constant depend on the " latitude " 
coordinate x (2> . Hence g* is a function of a: (1) alone whilst # 3 is 
a function of Xi and x$ alone. 

(c) ^4 does not involve the directional coordinates x (2) and a <3) 
and so is a function of z (1) alone. 

Accordingly, then, z (3) does not appear in the expression for 
(&)* and so, in addition to G u = 0, Gu = 0, G& = we have 
GIZ = 0, GK = 0. We must wiite down the five equations 
On = 0, On = 0, (? 22 = 0, G 3 3 = 0, G M = 0. The fact that HI 
is a function of z (1) alone and HI = 1 (a: (1) being the arc distance 
along the geodesic curves x (1) ) gives {14, 42} = and from 
Gu= {13,32} = Owe get 






which gives, on integration with respect to x (l \ 

= a function independent of 0:1 (A) 
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Go. = yields 



Qu = gives 



which on integration with respect to x 03 gives 



dH 

independent of a> (fi) 



Eliminating between (Q and U) wegetff s *^ i 



ent of x (1) . Since it cannot involve any variable but aP* 
have 

H^H* = a constant a, let us say; 

primes denoting differentiations with respect to x (1) . 
Again from (A) |^ - F 2 X a function of z< 2 > = ff 2 ^ say 

where < is a function of z (2) alone. Then jff 3 = fl"^ + / where 

1 $ s /7 

/ is a function of (1) alone. Now (5) shows that -- * is a 

1 



function of x (1) alone so that its derivative with respect to z (y > 
vanishes. Evaluating this derivative we find 

d<t> "-"H 2 } =0 



We can now proceed in various ways; either make 4> a constant 
or fH*' /'J? 2 (of which the second factor is the derivative) a 



constant giving / = const. X H^f -^-^ . We choose the latter 

/1 2 

alternative and make the constant zero so that / = giving 
# 3 = Hrf where <t> is a function of x (2) alone. 
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<p is determined by means of the equation <?ts 0, This gives 

f + *W-0 09 



On substituting jff $ = J3a# in (Z>) we find that - -^ is equal 



to a function of sf alone; but from its form and the definition 
of ^ it cannot involve a^ and so must be a constant. This 
coas^tant may, by a proper choice of unit for o; (2) , be put either 
1 or zero. We choose the first alternative and find, by suitably 
choosing the origin of measurement for o? (2) , <t> = sin a;< 2) . 
Gu = gives 

, H, 3*H* . aar, 



and on substituting 4 == sin s (2) , jff 3 = -^2^, both (D) and () 
yield the same equation 

M 2 " + H*H*' | g + ^ J = 1 (DO 

(B) gives 

TT rr TJ // 

2^- + ^- = (50 

Jfl2 124 

On differentiating ((?') and eliminating HzHt" we find 

jTI 2 "4 == -"2 "4 

which gives on integration 



where /5 is an arbitrary constant. 

Eliminating H 4 between (C f ) and (C 7 ') we have # 2 " = ^ 

f/2 

which on integration gives 



where 7 is an arbitrary constant. 
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Putting #4 = (ffsO/jS in (#') we have 



so that 1 = 2y giving 7 = f and hence finally ffj is determined 
by the differential equation 



and then 



sin 2 



It is usual to change the coordinate xP>, leaving the others 
unaltered. We write z (1) = z (1) (2/ (1) ) where y (1) = ff f . 



and we have 
- 



This is the form chosen by Einstein (that it is only one of many 
is evident from its derivation). If a& = it reduces to the 
well-known Euclidean form where y (l) = r, ar (2) = 8, # (3) = 4> are 
space polar coordinates. It remains to attach some physical 
significance to the constant a/3 and to take up the problem (6) 
stated at the beginning of this chapter. In order to conform 
to the usual notation we write henceforth y (l} = z'r; x ( ' 2) 8} 
z (3) = 0; z (4) = t where i 2 = - 1. 

Choosing the unit of x (4) or t so that /3 2 = + 1 and writing ia 
= a we have 

(&) 2 = - j ( 1 - ~ Y l (dr)^ + r 2 ^) 2 + r 2 sin 2 8(drf j 

+ (l - ? 
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3* DBTBSIGHATION OF THE PATH OF A FEEBLY MOVING 
PAETICLE 

A physical km c inertia is postulated to the effect that a 
/radjr moving material porttcfo in a gravitational field will follow 
the 3K>0-minimal geodesic lines of the four-dimermoTicd space time 
continuum which, for the single gravitating center, has the 
metrical gecsaetry characterized by the form given above by 

<*)*. 

A second postulate is that rays of light fottow the minimal 
P*xiwi--tlK>se for which ds = 0. In the ordinary Euclidean 
space these lines are imaginary, Le., have points with imaginary 
coordinates but the occurrence of the negative signs in the ex* 
presaion for (<&)* gives real minimal lines in our problem. For 
example, the light rays directed towards or away from the cen- 
ter, those for which B and <j> are constant, are characterized by 
the equation 



_ 1 _ =0 - 

In order, then, to solve the problem of the free motion of a 
material particle we have merely to determine the non-minimal 
geodesies whose equations are 

z<> + {lm, p}#*>df** = (cf. Ch. 6. 1) 

the dots denoting differentiations with respect to the arc length 
along the geodesic. For an orthogonal space of four dimensions 
these simplify to four equations of the type 

<+ {11, 1}< 2 + {22, l)(x< 2) ) 2 + {33, 1 }(*<>)* 

+ {44, l}(x>) 2 + 2{12, l}$<i< + 2{13, l}x<w 

+ 2{14, l)x (1) x (4 > = 

However we need use only three of these equations, replacing the 
fourth by r x (r) i ( * ) = 1 which is easily seen to be a consequence 
of the differential equations 

= 
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we multiply these by g n and use r as an mdbral symbol to 
;ain 





I then avail ourselves of the definition of the a?rabofc 
, ]. (Ch. 6, 1.) ^ (gwM.) is found to be sxo). In oar 



iblem it is convenient to omit the first of the 

: other three simplifying, on using the values for the three- 

ex symbols given (Ch. 7, 2), to 

where &.=*- \l--J ; 

*=-**; 

g* = r 1 sin* 9; 



these we have to add the first integral 

0ir 2 +0 2 2 +03? 2 +<M 2 = 1 (D) 

uations (B) and ((7) are immediately integrable giving 

<7s< = constant = h say (50 

1 

g4 = constant = + C say (CO 

on substituting the values of #3 and g 



r 2 sin 2 0= ft; l-yi= C 
uation (^4) may be written 

^ (r 2 ^) - r 2 sin 6 cos 0(() 2 = 

\ now proceed to eliminate the parameter s and find a relation 
meeting and 0. Assuming that < ^ (< = constant is a 
jcial case which is susceptible to the analysis given below on a 
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mere interchange of 8 and <p) we have 

so that 



+ W}<p - iV I? r'tf + 2 



Cte substituting the value of < from () we have 
- 

where we denote differentiations with respect to the new inde- 
pendent variable <}> by primes. Equating this to r 2 sin 6 cos 0($ 2 
and dividing out by r 2 ^ we obtain 

V 9 - 2 cot 0(0') 2 = sin 9 cos 
If now we choose our directional coordinate 6 so that initially 



We see that 9" = and then on differentiating the above equa- 
tion with respect to <p, 0'" = and so for all the other derivatives, 
i.e., 6 is a constant as <p varies. Otherwise expressed the general 
integral of the equation for as a function of <p is found by 
writing z = cot 6 yielding z" + z = to be cot 6 = L cos (< + M ) 
where L and M are arbitrary constants. We choose our initial 

conditions as above so that L = giving 6 = -r . Putting in this 

value for 9 we find 

A = h (B f ) 



(C") 

\ ' / 

and from (D) 
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Just as in the usual Newtonian treatment of planetary motion, 
it is convenient to write r = Iju and to again use ? g the 

independent variable. We have 



and then (ZX) yields, on making use of (C"), 



Now, in the Newtonian treatment, the equation giving the path 
of a particle under a central force is 



where F is the acceleration towards the center and 4 ss r* z 

4& 

is the constant of areas. Instead of this we have on differen- 
tiating the equation (E) just obtained 



u"+u ~ 

h 

so that we may, in a general manner, express Einstein's modifica- 
tion of the Newtonian law of gravitation by saying that there is 
superimposed to the inverse square law attraction an inverse 
fourth power attraction, the relative strength of the attracting 
masses being as 1 : 3A 2 . It remains to determine, at any rate 
approximately, the nature and magnitude of the constants a, h 
and C which arose in the integration of our differential equations. 
For large values of r, and therefore small values of u. the New- 
tonian law is a first approximation and so neglecting the term in 
if in the equation for u", a = F/zr = /i/w; /i being the gravita- 
tional constant and m the mass of the sun. Hence if we choose 
our unit of mass so that ju = 1 , a m, where now m is what is known 
as the gravitational mass of the attracting center (notice that we 
have identified, for small values of u, our r and 4> with the usual 
polar coordinates of Euclidean geometry). The velocity of light 

directly towards the attracting center is 1 r - and accordingly 
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our unit of time is such that for small values of u the velocity of 
Ug&t is unity; Le., if the unit of length be 1 cm., the unit of 
tfaac employed is l[c seconds where c = 3.10 10 . In the theory 
erf relativity there is no absolute distinction between space and 
time and so we refer to our time unit as one centimeter (1 cm. 
being tbe distance traversed by light in one time unit). It is 
to be observed that ii* Newtonian mechanics gravitational mass 
mhas dimensions Pf"* so that if L and T have the same dimen- 
sions a m has the dimensions of a length. The equation 

// i 7ft 

" + * = _ 

of the Newtonian theory yields 

u-~=Pcos(<-<o) 

where P and < are arbitrary constants of integration. 
Comparing this with the polar equation of a conic 
tu = 1 + e cos (I = semi-latus rectum, e = eccentricity) 

we have = / = A(l ) where A is the semi-major axis. 
If T is the period of revolution 



h = gXA^yf ellipse = 
whence 



where co is the angular velocity of the planet. This gives for the 
sun m = 1.47 kilometers or 1.47.10 5 cms. For the planets then 
w/r is a small quantity of the order 10~ 8 . In order to determine 
the constant C we differentiate 

?? 

u = ^ (1 + e cos <) 
and find 
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and comparing this with the equation () we ham 



It is to be observed that the values of m, C and k obtained fa 
this way are found from the Newtonian theory and ao are to be 
regarded as first approximations. In particular we have iden- 
tified the A of (') 
Accurately 

But 



so that neglecting quantities of the order ID"" 8 

^ _, ^ 
ck"^ "5" 

Substituting the expressions just obtained in () we have to 
integrate the first order differential equation 



m 2 (l 
d<f>J 

This equation defines u as an elliptic function of ^>; or inversely 
<p as an elliptic integral. It simplifies the algebra somewhat to 
write mu = x and to put m 2 /A 2 = a. We have already seen 

777 

that m 2 /^ 2 = - 5; so that if 6 is not very nearly equal to 
' 



unity a is a small quantity of the same order of magnitude as m A 
or 10~ 8 . Our equation is now 



- r) 

Now the discriminant of the literal cubic 

1 + aifl 2 + a& + a 3 = 
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is 



For the cubic cm the right-hand side of the equation giving (dvl&p}* 
this is 

4oV + &**(! - &*) ~ 108^(1 - ^) 2 . 



On account of the small magnitude of a this is positive, the first 
term being the dominant one. Hence the cubic has three real 
tools which we denote, in descending order of magnitude, by 
ti, **, ** When a * the roots are , 0, 0, and so we try first 
t ka and find i = (1 e) or (1 + 0) and then secondly 
SB |-f ia and find i = 2. Hence, to a first approximation, 
the three roots of the cubic giving (rftj/^) 2 are 9* = or(l <?); 
9g as a (i + ^); t>i = J 2a. Further since (cfe/d^) 2 cannot be 
negative in the problem D must lie between 3 and 2 or between 
t>i and + oo . As r does not tend to zero v does not tend to oo 
and hence t) lies between t>s and tJ2. We have 



r* dv 



The variable v oscillates between the values v& and t? 2 ; at these 
values dv/d<p = 0, so that has an extreme value; as v passes 
through the value v% retracing its values both dv and the radical 
change signs so that $ steadily increases. The change in <p 
between two successive extreme values of 2?, i.e., between peri- 
helion and aphelion of the planet, is 







V2(fl t?i)(tJ 



It is convenient to make a simple linear transformation of the 
variable of integration. Write v = a + bz and determine the 
coefficients a and 6 of the transformation so that to the roots 3 
and #2 of the cubic will correspond values and 1 of z respectively. 
The values are a = z? 3 ; 6 = 2 v* and then the third root Vi 
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goes over into 2 = ^ where * = **"** . The cubic 2( - 

( 02X0 03) transforms into 2Wa(l *) r 1 
so that 




This simplifies considerably on -writing z = sufj 




Now 

_ ^ ~ Pa _ 2oag 



ea 



(to a first approximation) is a small quantity of the same order 
of magnitude as a; hence we can expand (1 If sin* ffT 1 ** in a 
rapidly convergent series and a mere integration of the initial 
terms will give a very good approximation to k<p. The multiplier 
of the integral is 



2 ^ = 2/H - <>3 - 2[1 - 2a(3 - )]-'/ = 2[1 + a(3 - e)] 
and using J* * sin 2 6 d6 = ir/4 we find 



but fc 2 = 4a6 to a first approximation so 

A<p = TT{! + a(3 - e)}{l + ae} = 7rjl + 3a! 

Hence in a complete revolution the perihelion advances by an 
amount equal to 

O m 2 3m 

3a = 3 TT- = 



of a complete revolution, T being the period in our units. If we 
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wish to use the period in seconds and measure A in kilometers 
then the unit of time in the formula given is the time it takes 
light to trsvdi 1 kilometer 1/3.10* seconds; hence if T is the 
period in seconds the fractional advance of the perihelion per 

revolution is ^.^^^ - -z On substituting the values of 



A> f t aad * for Mercury's path this works out to be an advance 
of 43" per century. For the other planets e is much smaller than 
lor Mercury and the amount of advance of perihelion is much 
smaller; sare in the case of Mars the predicted advance is too 
small to be detected by observation. 

4. THE PATH OF A LIGHT RAT IN THE GEAVITATIONAL FIELD 
OP A SINGLE ATTBACTING CENTER 

These paths satisfy the equation (ds)* = or ds = 0; they 
are geodesies since, da being the non-negative root of the expres- 
sion for (&)*, no curve can have a negative length. The method 
of the preceding paragraph does not, however, immediately 
apply since the arc length a along a light ray, being a constant, 
cannot be used as an independent variable or parameter in terms 
of which the coordinates x may be expressed. Further in the 
discussion of Ch, 6, 1, it was assumed that the integral 



could be expanded in a Taylor series in powers of a. so that the 
existence of the derivative (d//da) as=0 was presupposed. It is 
apparent, however, on differentiation of 



ds = Tg lm dx ( dxW (/, m umbral) 

that if ds = when a = 0, (ds) becomes meaningless when 

da 

a = on account of the zero factor (ds) ASs0 which occurs in the 
denominator. These difficulties are overcome in the following 
manner. If we investigate those curves (non-minimal) for which 
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the first variation of the integral I s= f(&f* & oro we are kd 
to exactly the same differential equatiom as theee of QL &, { 1, 
which express the fact that the first variation of I J*4* fa aoro! 
Accordingly we now derive the equations of the mmm^ geodcstea 
from the fact that the first variation of I >*(&)* is **o, A 
being zero along the curves. The coordinates 2 arc imyyjjstJ 
expressed in terms of any convenient parameter aad differentia 
tions with respect to this parameter are denoted fey 
The Euler-Lagrangian equations are (ef. du 6, 1) 

dF 



where 

Fm(d,ym ^rfW'asW* ft umbl) 

Hence 



or 

ffrix (l)// + [lm, r]*> V->' = ft m umbml) 

Multiplying by ^^ and using r as an umbral symbol we obtain 
x ( &" + [lm, p}x^'x (m)/ = (p = 1, 2, 3, 4) 

which are exactly the equations of Ch. 6, 1. The first integral 
of these equations which has already been mentioned may be 
very briefly obtained as follows. Since F s gi m x {I) 'x (m >' is 

* Attention should, however, be called to the fact that this integral is not, 
properly speaking, a line integral at all; its value depends not only on the 
curve over which it is extended but on the particular parametric mode of 
representation chosen for this curve. In order that the value of the integral 
should not be dependent on the parametric representation the integrand should 
be positively homogeneous of degree unity in the derivatives x'. 
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of degree 2 in the x 7 we have, by Euler's theorem 



on homogeneous functions, 



Ik result imniedbtely verifiable directly (r umbral). On multiply- 
ing the equations 

9F _ d ( dF \ Q 



by a0*' and using r as an umbral symbol we obtain 

*.-*.(&>*. 

d*\ da:"' 



or 

dv dv 

showing that F is constant along the geodesies. The constant 
is now zero instead of unity as it was in the case of the non- 
minimal geodesies. 

Before proceeding to calculate the deflection of the light rays 
it will be well to prove an often quoted property of them. In a 
statical gravitational field the time coordinate x (4) does not enter 
into F s (cfe) 2 . Hence 



dF \ 

- , 1 

dx^ ) 



A 

or - = const. 



If, now, in the discussion of Ch. 6, 1, instead of keeping both 
ends of the " varied curve " C(or) fixed, we had allowed the ends 
to vary also, the part of SI which came outside the sign of integra- 
tion when we integrated by parts would not vanish automatically. 
Since the first variation is to vanish when the end points are 
fixed as well as when they vary the part under the sign of integra- 
tion vanishes as before yielding the Eulerian equations but in 



PBDBLEMS IN 
addition we have the end condition 



If now all the coordinates but z^ are fcept fixed 

&e<i> = o = &c^ fa^ 

aip 

and we find since - ,is constant over tibe extremal 



and as 



we have 

= 



which is known as the Fermat or Huyghtns 1 Principle of Lecut 
Time. It is an immediate consequence of the absence of y?& 
from (<fo) 2 ; there is a similar theorem for the symmetrical 
attracting center: 





but this has no special utility. The Fermat Principle states that, 
given two fixed points in space (by fixed is meant that the three 
space coordinates for an observer attached to the gravitating 
center are constant), a light signal passes from one to the other 
in such a way that the first variation of the time interval is a 
minimum. 
With the same notation as that employed in 3 we find 

9& f = - A; gjf = C 

where h and C are constants and we find exactly as before that a 
proper choice of our initial conditions for 6 enables us to 
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write 9 m */2. The only difference is that (DO is replaced by 

(1 - 2m/rrW + fW - (1 - 2/r)(O f - 
whence on writing r = 1/t* and using rV = * we find 



In order to get an idea of the order of magnitude of the constants 
C, A of integration we make a first trial-approximation. The 
largest value that t* can have is IfR where R, the radius of the 
sun, 897,000, the units being kilometers. Hence we neglect, 
for tiie moment, the if term in comparison with the others and 
find at once 

Q 

u = 7- sin (tp <po) 
fl 

where ^ is a constant of integration. Hence C/h is the largest 
value of u and is therefore a small quantity of the order 1/10 5 . 
Denoting this small quantity by a (a is the positive square root 
of C*/h*) we have 



The discriminant (cf. Ch. 7, 3) of the cubic on the right is 
4o 2 (l 27m 2 o 2 ), a positive quantity, so that the three roots are 
real. When a = they reduce to l/2m, 0, 0, so that trying in 
turn ka and (1 /2m) + ka we find the first approximation to the 
three roots u^ = or, uz = a, u\ = l/2m where we have ar- 
ranged the roots so that u^< u^< u\. For a second approxi- 
mation, we try in turn a + ko?, a + Aro 2 , (l/2m) + Aro 2 for 
^3, H$, "1 respectively and find 

u* = a + ma 2 ; Wo = a + ma 2 ; ^! = (l/2m) 2ma 2 . 

We now, as before, determine a linear transformation which 
sends u = tx 8 into z = 0; w = i^ into 2=1. It is u = a + bz 
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iere a = wa> b = u^ u z and then the third toot * i 
es into z = 1/& 2 where 



ow the cubic 2?nzj 3 w 2 + o? cannot be negative in oar problem 
>r can u itself. At remote distances from the son * no 
At initially u = and it increases to u = % at which point * 

is a maximum value, since ( -~ J = there. Then *i begins to 
^crease and the radical *>l2mtf u 2 + a? in the expression for 



, 
ay 



so changes sign so that d<p keeps its sign. The angle between 
point at a remote distance and the perihelion of the light ray 

given by the integral | The excess of 

Jo V2mw 3 if + of 

vice this over TT is the deflection D experienced by the ray. 
ence 



^mu - u* + of 
hich on writing ?/ == a + 62 becomes 



2) f rj 2 j 



n making the final substitution 2 = sin 2 6 this becomes 
D + ,r 



^P Vl-A-sin 2 ^ 

'ow 

= 4ma + higher powers in a 



~ + a - 
2m 
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so that, m being 1.47, if is a small quantity of the same order 
as ou Hence (I W sin* fl)" 4 ' 1 can be expanded in a rapidly 
convergent series and an integration of the initial terms of this 
series gives a high approximation to D + T. On substituting 
the values of k and b the multiplier of the integral becomes 



4(1 -f 2ma)- 1/1 - 4(1 - ma) 



whilst the lower Hmh of the integral is 



^^-sin-^d-ima) 




Here it is necessary to use the second approximation since u* is 
to be divided by u t * itself a small quantity of the first order. 
On expanding 



by Taylor's theorem we have for the lower limit (:r/4) f ma 
so that 

D + T 4(1 -wajftf + Y (0- sine cos 0) + 



In the term multiplied by 



it is sufficient to take the rough approximation ?r/4 to the lower 
limit and we have 



gvng 

D 

a, in this expression, is the maximum value u^ of u (to a first 
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approximation), Le., is the reciprocal of the radius of the * 
An idea as to the closeness of this approximation is obtained 
by using the second approximation 



mc? 
The positive root a of this quadratic is 



so that writing a = l/R is equivalent to neglecting mfR in 
comparison with unity or to a neglect of 1 part in 5.1(P. OB 
substituting m = 1.47, R = 697,000 in the expression D 4m/B 
and converting this radian measure into seconds of arc we find 
the value 1.73" predicted by Einstein for a light ray which just 
grazes the sun.* 

*For a fuller discussion of the problems dealt with in this chapter refereooe 
is made to two papers by the author in the Phil. Mag. of dates Jan. (1922) 
and March (1922) respectively. For an alternative treatment of the subject 
matter of 2 the reader should consult the paper Concomitants of Quadratic 
Differential Forms by A. R. Forsyth in the Proc. Roy. Soc. Edin. May (1922) . 
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